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ABSTRACT 


A  methodology  is  described  for  identifying  least-cost  forest 
harvesting  systems  using  shortest  path  network  analysis.  Components 
of  the  harvesting  system  are  depicted  in  the  form  of  a  directed,  acyclic 
network  describing  the  sequence  of  alternative  methods  for  satisfying 
the  system's  objectives.  Nodes  in  the  network  represent  completed 
operations  (events).  Arcs  joining  the  nodes  have  values  (durations) 
representing  the  cost  of  the  event.  The  optimal  (least-cost)  and 
near- optimal  forest  harvesting  systems  are  found  by  selecting  and 
ranking  by  lowest  total  cost,  different  paths  through  the  network. 

A  package  of  computer  programs,  written  in  BASIC  for  a 
Hewlett-Packard  9845B  mini-computer,  was  developed  for  implementing 
the  methodology.  Routines  aid  in  the  definition  of  the  network,  and 
input  of  variables  for  each  activity.  Activity  durations,  measured 
in  dollars  per  cubic  metre  of  wood  produced,  were  calculated  from 
estimates  of  hourly  costs  and  production  for  specified  conditions. 

A  before- tax  cash  flow  model  calculates  activity  costs  on  an  annual 
equivalent  cost  basis.  Productivity  is  determined  using  prediction 
equations  to  estimate  activity  cycle  times  for  the  conditions. 

Sample  application  of  the  methodology  is  made  to  identify 
the  cheapest  system  for  logging  and  transporting  tree-length  material 
to  a  processing  plant  under  Eastern  Canadian  conditions.  The  optimal 
system,  a  swather-type  feller-primary  transporter  presently  under 
development  combined  with  roadside  flail  limbing,  was  significantly 
cheaper  than  alternative  systems.  The  optimal  system  selected  was 
found  to  be  insensitive  to  fuel  prices  with  fuel  prices  up  to  double 
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present  costs.  Sensitivity  analysis  was  also  used  to  identify  the 
least-cost  harvesting  systems  under  a  range  of  tree  diameters  and  stand 
volumes.  The  swather  system  was  best  under  most  tested  conditions. 
However,  with  large  tree  sizes  and  low  stand  volumes,  a  system  with  a 
tractor-mounted  shear,  tracked  grapple  skidding  and  roadside  flail 
limbing  produced  the  lowest  costs.  A  manual  felling  system  produced 
the  lowest  costs  with  medium  sized  trees  (24  cm  DBH)  and  low  stand 
volumes. 

Shortest  path  network  analysis  was  found  to  be  a  useful 
method  of  designing  forest  harvesting  systems.  The  ability  to 
identify  both  optimal  and  near-optimal  systems  is  important  when  the 
differences  between  most  systems  are  as  small  as  those  observed  during 
the  case  study.  The  technique's  biggest  advantage  is  that  many  more 
alternatives  than  would  usually  be  evaluated,  can  be  examined  quickly 
and  in  a  manner  which  permits  easy  comprehension  by  non-technical 
personnel. 
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1.  INTRODUCTION 


Since  the  invention  of  the  axe,  increasingly  sophisticated 
machines  and  systems  have  been  developed  for  harvesting  wood  fibre 
from  forested  lands.  Logging  management  personnel  must  select  that 
combination  of  harvesting  equipment  which  is  best  suited  to  the 
particular  operating  environment  from  the  array  of  alternative 
harvesting  systems  which  have  been  or  could  be  developed1. 

The  best  system  for  any  specific  operating  environment:, 
whether  physical,  economic  or  social,  can  only  be  identified  after 
all  possible  methods  of  meeting  the  objectives  of  that  system  have 
been  investigated.  The  research  discussed  herein  investigated  the 
use  of  a  computer  program  which  selects  and  ranks  alternative 
combinations  of  machines,  facilitating  the.  selection  of  equipment 
and  the  design  of  the  least-cost  forest  harvesting  system. 

1.1  Research  Objectives 

The  objectives  of  the  research  undertaken  for  this  thesis 
were  as  follow: 

a)  To  describe  with  directed  networks  possible  alternative 
forest  harvesting  systems; 

b)  To  compile  engineering,  economic  and  time  study  data  on 
selected  logging  equipment; 


1  It  is  suggested  that  readers  unfamiliar  with  forest  harvesting  terms 
consult  a  standard  reference  on  terminology  e . g .  L108;  144;  145]. 
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c)  To  develop  production  cost  prediction  equations  for 
harvesting  equipment; 

d)  To  demonstrate  the  use  of  shortest  path  network  analysis 
as  a  means  of  selecting  optimal  forest  harvesting  systems 
and 

e)  To  evaluate  the  sensitivity  of  the  logging  system  to 
variations  in  its  operating  environment,  and  so  determine 


optimum  operating  zones. 
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1.2  Environment  Facing  Canadian  Loggers 

Nadeau  [111],  Morrison  [110],  and  others  [68;  123;  142]  have 
shown  a  deterioration  in  the  competitiveness  of  the  Canadian  forest 
products  industry  compared  with  most  other  producing  countries. 

Seventy  to  eighty  percent  of  the  differences  in  total  pulp  and  paper 
manufacturing  costs  between  Canadian  and  American  producers  has  been 
attributed  to  high  wood  costs  [142,  p.  41].  Although  the  situation 
is  less  critical  in  the  lumber  industry,  some  regions  of  Canada  have 
competitive  disadvantages  for  the  same  reason.  These  high  wood  costs 
are  attributable  to  a  number  of  problems  within  the  physical,  economic 
and  social  environments  facing  Canadian  loggers. 

Increasing  demand  for  forest  products  has  resulted  in  greater 
harvests  as  shown  in  Figure  1,  and  forced  greater  utilization  of  what 
had  earlier  been  considered  stands  and  species  of  marginal  value. 
However,  Reed’s  [138]  study  of  Canada’s  reserve  timber  supply  concludes 
that  the  economically  exploitable  forests  are  now  almost  completely 
allocated,  and  only  improved  utilization  and  intensified  forest 
management  will  permit  further  expansion  of  the  industry. 

Reduced  average  tree  size  and  rapidly  rising  wages  have 
prompted  a  trend  towards  the  complete  mechanization  of  forest 
harvesting  operations  as  a  means  of  raising  productivity.  Hand 
saws  and  horsedrawn  sleighs  have,  for  the  most  part,  been  replaced 
by  power  saws,  skidders  and  other  mechanized  equipment. 
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Figure  1.  Total  forest  production  in  Canada  1960  -  1978  [147;  149] 


The  urbanization  of  society  reduced  the  labour  supply  of 


an  industry  which  had  been  labour  intensive  prior  to  1960.  The 
increased  labour  productivity  offered  by  mechanization  became  a 
necessity  under  these  circumstances.  Table  1  shows  that  the  average 
annual  productivity  of  forest  workers  doubled  in  the  period  from 
1960  to  1978. 

Labour  costs  also  increased  during  the  period,  both  as  a 
result  of  inflation  and  the  need  for  a  wage  and  fringe  benefit  package 
to  neutralize  the  urbanization  trend.  In  terms  of  constant  value  1971 
dollars2,  average  wages  more  than  doubled  from  $4,830  to  $10,039  between 
1960  and  1978. 

Over  the  shorter  term,  the  strong  inflationary  pressures 
following  the  energy  crisis  in  1971  produced  dramatic  increases  in  the 
costs  of  energy,  lubricants  and  purchased  supplies  as  Figure  2  shows. 

These  factors  combined  to  produce  the  equivalent  of  a  325% 
increase  in  the  total  cost  per  cubic  metre  on  a  current  dollar  basis 
between  1960  and  1978  as  shown  in  Figure  3.  On  an  absolute  basis, 
average  unit  costs  rose  at  an  average  annual  rate  of  over  2%  during 
the  period. 


2 


Constant  value  1971  dollars  were  used  to  remove  the  inflationary 
effects  during  the  period.  Various  price  indices  have  been 
applied  to  some  cost  figures  as  outlined  in  Appendix  A. 
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Table  1.  Annual  Size,  Productivity  and  Average  Wage  for  the  Labour  Force  Employed  in 
Canadian  Logging  1960  -  1978  [147;  149]. 
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Figure  2.  Breakdown  of  total  logging  costs  in  Canada  1960  -  1978 
[147;  149;  150]. 
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per  cubic  meter  in  Canada  1960  -  1978  [147;  149;  150]. 


Figure  3. 


Although  mechanization  has  solved  some  of  the  problems 


associated  with  labour,  it  has  created  many  of  its  own.  The. 
complexity  of  logging  machines  has  increased  the  necessity  for  the 
training  of  woodlands  employees.  Higher  demands  for  skilled  labour 
has  often  placed  the  industry  in  direct  competition  with  urban 
employers  for  trained  personnel.  With  mechanization,  machine 
maintenance  has  become  an  increasingly  important  area  of  concern. 
Organizational  restructuring  and  the  development  of  new  infra-structure 
became  necessary  to  maintain  and  control  logging  machines. 

Moreover,  the  sheer  magnitude  of  the  multi-million  dollar 
investments  in  the  capital  assets  required  for  large-scale  harvesting 
operations  has  become  a  major  problem.  The  size  of  the  actual 
investment  varies  with  the  degree  of  mechanization  and  the  amount  of 
wood  cut,  but  as  Figure  4  shows,  capital  expenditures  have  risen. 

The  ratio  of  capital-to- labour  input  costs  increased  by  20%  in 
constant  dollars  between  1960  and  1978.  This  suggests  that  the 
doubling  of  man-year  productivity  shown  in  Table  1  is  the  result  of 
a  substitution  of  capital  for  labour  during  the  period. 

To  measure  the  relative  effectiveness  of  these  production 
factors,  a  comparison  of  the  productivity  indices  for  capital  and 
labour  expenditures 3  is  made  in  Table  2.  During  the  time  frame, 
capital  productivity  decreased  slightly,  while  there  was  a 
compensating  marginal  increase  in  labour  productivity. 


The  productivity  indices  were  calculated  using  the  equation 


Index  =  — 


Production  Implicit  Price  Index 


x 


Expenditures  “  100 

using  expenditures  on  capital,  labour  and  a  combination  of  both. 


. 


»- 


. 
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Figure  4.  Total  capital  invested  by  the  Canadian  forest  products 
industry  in  logging  construction,  equipment  and  repairs 
between  1960  -  1978  [150]. 
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However  since  capital  was  being  substituted  for  labour, 
the  combination  of  these  two  factors  gives  a  more  accurate  estimate 
of  productivity.  This  approach  indicates  that  a  0.1%  annual  decrease 
in  the  overall  productivity  of  the  Canadian  logging  industry  occurred 
between  1960  and  1973.  While  new  technologies  were  twice  introduced 
into  the  industry  during  these  years  (skidders  and  harvesting  equipment) , 
they  had  no  major  impact  on  over-all  productivity.  This  supports 
Morrison's  conclusion  that  "neither  capital  nor  labour  were  used 
effectively"  in  the  transition  from  a  labour  intensive  to  a.  capital 
intensive  industry  [110,  p.  vi]. 

The  net  result  of  this  stable  productivity  has  been  rises 
in  wood  costs  to  the  point  where  United  States  producers  generally 
enjoy  a  lower  material  costs  than  all  Canadian  production,  regions 
other  than  Interior  British  Columbia  [123,  p.  11-26], 

This  has  been  a  main  contributing  factor  to  low  the 
profitability  that  threatens  the  availability  of  capital  in  the 
future.  Since  the  threat  of  shortages,  small  free  sizes  and  relative 
inaccessibility  compared  to  competitors  will  continue  if  not  worsen 
in  the  future,  one  of  the  main  problems  now  confronting  the  Canadian 
logging  industry  involves  improving  the  productivity  of  capital  and 
labour  expenditures.  This  research  investigates  a  method  which  aids 
woodlands  management  personnel  in  the  analysis  and  select  of  systems 
most  likely  to  improve  productivity. 

A  brief  outline  of  the  basic  forestry  harvesting  systems 
and  their  component  subsystems  follows. 


' 


1.3  Forest  Harvesting  Systems 

Nadler's  general  definition  of  any  system  [112,  p.  41]  as 
the  specified  and  organized  conditions  for  the  elements  of  function, 
inputs,  outputs,  sequence,  environment,  physical  catalyst,  and  human 
agents  detailed  for  each  element  in  physical,  rate,  control,  and  state 
dimensions  clearly  applies  to  harvesting  systems. 

These  elements  relate  to  forest  harvesting  systems  as  follows. 
The  function  is  to  harvest  and  transport  a  quantity  of  wood  fibre  to 
some  location  for  additional  processing  at.  minimal  cost.  Large 
quantities  of  trees  and  energy  are  the  system's  main  inputs.  The 
output  of  a  logging  operation  is  wood  fibre  at  the  consuming  mill's 
site.  Between  the  system's  inputs  and  outputs,  the  operation  follows 
a  sequence  of  processing  steps.  The  conditions  within  the  forest, 
society  and  the  economy  are  the  environment  in  which  the  system 
operates.  The  physical  catalysts,  which  aid  the  conversion  of  inputs 
to  outputs,  are  the  various  logging  machines.  Woodlands  workers  and 
management  are  the  human  agents  associated  with  the  system.  Thus 
forest  harvesting  operations  are  akin  to  other  systems.  As  such, 
systems  analysis  is  one  method  of  analyzing  logging  systems  and 
identifying  the  optimal  system. 

There  are  five  distinct  general  categories  of  forest 
harvesting  systems: 

1.3.1  multiple-length, 

1.3.2  tree-length, 

1.3.3  full-tree, 

1.3.4  whole- tree,  and 


1.3.5  chipwood. 


> 


I  & 


■ 


■ 


. 


Figure  5  illustrates  the  historical  use  of  the  different 
systems  in  Eastern  Canada.  Note  that  the  basis  of  comparison  is 
wood  form  upon  delivery  to  roadside 5  rather  than  form  upon  delivery 
to  the  mill  as  used  elsewhere  in  this  thesis. 


MULTIPLE  LENGTH 

TREE  LENGTH 

FULL  TREE 


Figure  5.  Historical  use  of  different  forest  harvesting  systems 
in  Eastern  Canada  [104,  p.  5]. 
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1.3.1  Multiple-Length  Harvesting  Systems 


Multiple-length  harvesting  systems  involve  the  production 
and  transport  of  logs  of  less  than  the  merchantable  height  of  the  tree 
to  the  processing  mill.  Generally  if  the  logs  are  over  three  meters 
in  length,  the  system  is  called  longwood,  and  if  under  three  meters, 
shortwood.  Both  types  involve  the  same  basic  operations  but  with  some 
differences  due  to  the  log  length.  The  main  subsystems  are  as  follows 

a)  felling  of  the  tree, 

b)  bolt  preparation  i . e .  limbing,  topping,  and  bucking, 

c)  transportation  of  wood  from  the  stump  to  a  central  point 

d)  loading, 

e)  transportation  to  the  processing  mill,  and 

f)  unloading. 

The  sequence  of  the  subsystems  is  variable,  with  the  possibility 
of  the  transportation  from  stump  to  landing  and  bolt  preparation 
subsystems  being  interchanged.  It  is  also  common  to  have  bolt 
preparation  disunited,  with  for  example  limbing  and  topping  done 
at  the  stump,  and  bucking  at  the  landing. 

Historically  in  Canada,  the  multiple-length  system  has  been 
the  primary  method  of  harvesting  wood  fibre.  The  ease  of  handling 
with  short  material  when  technology  was  less  developed,  the  require¬ 
ments  of  groundwood  pulp  mill  equipment,  the  delivery  to  the  mill  of 
only  "acceptable"  portions  of  the  tree,  and  legislated  limits  on  the 
widths  and  lengths  of  truck  loads  were  the  major  justifications  for 
the  continued  usage  of  the  system.  Technological  advances  in  the 
last  25  years  have  largely  negated  these  advantages. 
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1.3.2  Tree-Length  Harvesting  Systems 

A  tree-length  system  produces  and  transports  wood  in  logs 
equal  in  length  to  the  merchantable  height  of  the  tree.  The  actual 
length  varies  with  the  forest  environment,  tree  species’  characteristics 
and  the  merchantability  limits,  but  the  transported  logs  usually  exceed 
12  meters  in  length  and  rarely  are  longer  than  18  meters.  This  system 
involves  the  same  subsystems  as  the  multiple-length  system,  but  the 
bolt  preparation  stage  consists  only  of  limbing  and  topping. 

1.3.3  Full-Tree  Harvesting  Systems 

Full-tree  logging  results  in  the  delivery  of  a  felled  tree 
with  the  limbs  largely  intact  to  the  processing  plant.  The  bolt 
preparation  subsystem  is  completely  removed. 

1.3.4  Whole-Tree  Harvesting  Systems 

The  whole  or  complete-tree  harvesting  system  involves  the 
harvesting  of  all  wood  fibre  within  the  tree,  from  roots  to  branches. 
There  are  potential  problems  with  nutrient  cycling,  soil  erosion  and 
regeneration  associated  with  the  system,  so  it  remains  largely  unused 
in  Canada.  However,  Keays  [71;  72;  73;  74;  75]  suggested  that  the 
system  has  great  potential  for  some  products,  since  the  utilization 
of  the  total  wood  fibre  produced  by  the  tree  could  be  increased  by 
as  much  as  forty  percent  for  most  Canadian  species. 

1.3.5  Chipwood  Harvesting  Systems 

Chipwood  harvesting  systems  deliver  wood  chips  to  the 
consuming  mill.  The  chipping  of  trees  in  various  forms  is  possible, 
but  the  inclusion  of  tree  bark  with  the  chips  creates  problems  with 
the  pulp  quality  with  present  pulping  processes.  The  system’s  output 


are  also  unsuitable  for  utilization  in  lumber  and  veneer  mills ,  so  the 


■ 


> 
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potential  usage  of  the  system  is  limited  to  pulp  and  particle  board 
products.  The  purest  form  of  this  system  would  be  chipping  of  the 
whole-tree  at  the  stump,  and  subsequently  transporting  wood  chips  to 
the  mill. 

Having  described  the  basic  forest  harvesting  systems,  the 
literature  covering  different  techniques  which  have  been  applied  to 
the  analysis  of  these  systems  will  be  reviewed  in  the  next  section. 
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2.  LITERATURE  REVIEW 

In  the  past  three  decades,  the  trend  in  forest  harvesting 
research  has  been  away  from  the  study  of  single  machines  to  the  study 
of  groups  of  machines  combined  into  systems.  Prior  to  the  mid-1960's, 
most  studies  investigated  harvesting  as  single  entities,  although  it 
was  often  recognised  that  many  situations  involved  several  activities. 
Since  that  time,  the  systems  approach,  studying  problems  as  a  collection 
of  activities  which  interact,  has  become  increasingly  common. 

Although  many  analytical  techniques  have  been  used  to  analyze 
harvesting  machines  and  systems,  the  application  of  a  particular 
technique  has  often  not  been  limited  to  one  use.  For  example,  a 
simulation  model  could  be  developed  for  designing  a  system,  and  could 
also  be  used  for  planning  and  control  purposes  once  the  system  existed. 
As  such,  the  techniques  described  in  the  literature  have 
been  separated  into  their  primary  area  of  usage  as  follows : 

2.1  Determining  Machine  Productivity 

2.2  Planning  and  Controlling  Logging  Operations 

2.3  Designing  Logging  Machines  and  Systems 


2.4  Synthesis 


. 


. 
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2.1  Determining  Machine  Productivity 

Most  analyses  of  logging  machines  and  systems  used  time 
study  data,  alone  or  in  conjunction  with  various  measurements  of  the 
operating  environment,  to  determine  machine  productivity  and  cycle 
times.  However,  data  analysis  has  become  increasingly  complex. 

Simple  statistics  such  as  means  and  standard  deviations  have  been 
replaced  by  multiple  regression  and  time  series  analysis.  In  recent 
years,  operations  research  models  using  probability  distributions  for 
time  elements  have  become  common. 

The  main  objective  of  most  early  studies  was  the  determina¬ 
tion  of  productivity  and  the  influence  of  environmental  factors  on  it. 
Research  such  as  that  by  McCraw  [100],  Bennett  et  al .  [8],  Dibblee  [35] 
and  Harvey  [63],  concentrated  on  detailed  time  studies  of  wheeled 
skidders.  F.esults  of  these  studies  Indicated  that  about  60%  of  the 
variation  in  productivity  could  be  explained  by  easily  measured 
variables  such  as  load  volume  and  skidding  distance.  Several  authors 
e. g.  [101;  9;  27]  attempted  to  improve  the  accuracy  of  predictions  by 
including  such  factors  as  surface  roughness,  slope,  brush  density  and 
ground  bearing  strength  but  with  little  success.  Variable  skidding 
crew  motivation,  changes  in  their  physical  effort  (rating)  in  response 
to  work  conditions,  and  interaction  between  all  these  factors  caused 
measurement  problems. 

More  recent  studies  of  skidding  operations  have  been  unable 
to  improve  the  accuracy  of  these  early  prediction  equations.  Studies 
as  recent  as  1976  by  Matthes  et  al .  [98]  explained  about  65%  of  the 
productivity  variation  using  a  six  variable  multiple  regression 
equation.  In  comparison,  McCraw’ s  and  Hallet’s  studies  [101]  made 
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during  the  1960's,  used  three  variables  to  explain  55%  of  the  variation. 
Cottell  et  al.  [27;  28]  had  accuracies  of  80%  and  more  using  equations 
with  two  independent  variables  and  shift-level  data. 

Kroger's  study  [77]  is  an  exception  to  these  relatively  low 
accuracies.  He  reports  coefficients  of  multiple  determination  exceeding 
94%.  However,  his  models  are  quite  complex  and  their  data  requirements 
make  them  of  questionable  value  for  practical  uses.  The  trade-off 
between  accuracy  and  usefulness  has  led  many  researchers  to  avoid 
complex  productivity  prediction  equations  with  higher  accuracy  since 
the  numerical  skills  among  those  using  the  equations  were  frequently 
poor  and  problems  with  measuring  factors  such  as  brush  density  and 
surface  roughness  consistently. 

In  more  recent  years,  both  detailed  and  gross  time-study 
procedures  have  been  used  for  logging  machine  studies.  The  machine 
evaluation  studies  made  by  Skogsarbeten  (Swedish  Logging  Research 
Foundation),  the  American  Pulpwood  Association,  the  Pulp  and  Paper 
Research  Institute  of  Canada  (PAPRICAN) ,  and  the  Forest  Engineering 
Research  Institute  of  Canada  (FERIC)  use  short-term  (3  to  5  day) 
detailed  studies  to  prepare  productivity  estimates  [19;  128;  84;  46]. 
Other  organizations  around  the  world  use  similar  techniques  for  the 
same  purpose,  e. g.  [10;  12;  60]. 

Folkema  [47]  and  others  at  FERIC  have  used  gross  time  study 
procedures  in  conjuction  with  time  recording  instruments,  i . e .  Servis 
recorders,  to  determine  long-term  machine  time  data,  productivity  and 


causes  of  machine  downtime. 


. 
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The  failure  of  the  early  researchers  to  quantify  the 
unexplained  variation  led  to  the  development  of  gross-data  procedures 
which  used  shift-level  or  day-long  periods  as  the  basis  of  prediction. 
Productivity  estimates  with  sufficient  accuracy  for  production  control, 
planning  and  budgeting,  and  in  some  cases  greater  accuracy  than  detailed 
studies,  were  found  to  be  obtainable  using  these  gross  data  procedures 
[28;  158].  Detailed  time-element  studies  have  been  reserved  for  method 
studies  aimed  at  improving  work  techniques,  the  evaluation  of  new 
machines,  and  for  determining  the  probability  distributions  and 
regression  equations  used  in  simulation  models. 

Two  conclusions  are  apparent  from  the  studies  chiefly  concerned 
with  determining  the  productivity  of  logging  machines  and  systems: 
machine  and  system  performance  has  a  stochastic  nature,  and  the 
interactions  between  system  components  were  usually  not  studied. 

2. 2  Planning  and  Controlling  Logging  Operations 

The  results  of  many  early  time  studies  were  intended  for 
use  in  the  planning  and  control  of  logging  operations.  Cunia  [31] 
and  Lussier  [91]  used  statistical  data  gathered  from  time  studies 
of  equipment  and  crew  performance  to  prepare  production  standards 
and  control  charts,  and.  as  a  basis  for  budgeting  and  the  evaluation 


of  future  systems. 
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2.2.1  Nomograms 

Many  authors  have  avoided  the  problem  of  low  numeracy 
among  potential  users  by  providing  nomograms  for  productivity  and 
cost  calculations,  e. g .  Ager  [1],  and  Legault  [85].  The  productivity 
curves  are  usually  established  using  regression  equations  developed 
from  time  studies.  Costing  models  of  various  types  have  been  used 
for  determining  operating  cost  of  the  equipment. 

The  so-called  "Weak  Link"  analyzer  developed  by  Baumgras  and 
Martin  [6]  uses  nomograms  to  aid  loggers  in  analyzing  their  harvesting 
systems.  The  productivity  and  costing  nomograms  presented  for  the 
felling,  skidding  and  transportation  subsystems  represent  a  complete 
systems  approach.  The  production  nomograms  are  limited  to  conditions 
in  the  Appalachian  region  of  the  United  States,  but  the  approach  is 
adaptable. 

The  economic  analysis  of  20  logging  systems  used  in  Quebec 
made  by  Conseillers  en  gestion  des  forets  (COGEF)  contains  nomograms 
for  adjusting  the  observed  productivity  of  logging  machines  [55]. 
Logging  chance  was  classed  by  degree  of  difficulty  1 . e .  easy,  medium 
and  inexploi table,  and  combined  with  operator  and  task  characteristics 
to  create  an  adjustment  factor.  Although  intuitively  the  curves 
appear  reasonable,  quantitative  support  is  absent. 
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2.2.2  Differential  Calculus 

Probably  the  earliest  application  of  a  comprehensive 
mathematical  analysis  to  the  planning  of  timber  harvesting  operations 
was  Matthew’s  treatment  of  road  and  landing  spacing  [99].  His  approach 
used  the  method  of  equal  areas  to  derive  the  average  skidding  distance 
for  regularly  shaped  areas  like  circles,  triangles  and  rectangles. 

The  average  skidding  distance  was  used  in  a  break-even  analysis  which 
equated  variable  skidding  costs  and  road  construction  costs  to  determine 
optimal  road  and  landing  spacings. 

Several  other  approaches  have  been  developed  for  investigating 
the  optimal  spacing  problem.  Lussier  [92]  revised  Matthew’s  work, 
recognizing  the  influence  of  terrain  conditions  on  skidding  productivity. 
Suddarth  and  Herrick  [153]  developed  formulae  based  on  integral  calculus 
for  the  average  skid  distance.  Their  work  demonstrated  that  Matthew’s 
average  skid  distance  formula  for  a  rectangle  was  incorrect.  Lysons 
and  Mann  [93]  showd  that  Matthews’  formula  for  circle  wedge  shapes 
was  also  in  error. 

All  these  approaches  assume  that  the  cutting  block  has  a 
regular  shape,  two-way  skidding  to  the  landing  is  possible,  variable 
skidding  costs  are  equal  from  both  sides,  and  that  stand  density  is 
constant  throughout  the  cut  block.  In  practice,  none  of  these 
assumptions  is  valid.  Additionally,  topographic  features  and  soil 
conditions,  such  as  streams,  muskegs  and  cliffs,  may  prohibit  positioning 
the  landing  in  the  central  location,  and  increase  the  actual  skidding 
operation,  rather  than  reduce  productivity. 


. 
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Peters  and  Burke  [125]  extended  the  integral  calculus 
approach  to  irregularly  shaped  areas,  using  a  procedure  requiring 
a  digitizer.  Peters  [126]  also  formulated  a  generalized,  direct 
solution  method  which  uses  Matthew's  expression  for  total  harvesting 
cost  and  the  integral  form  for  average  skidding  distance.  The 
results  showed  that  a  necessary  condition  for  minimum  costs  is 
that  variable  skidding  costs  equal  road  costs  plus  twice  landing 
costs . 

Donnelly  [37]  developed  a.  program  for  a  hand-held  program¬ 
mable  calculator  for  calculating  average  skidding  distance  for  any 
shape  of  cut-block.  His  method  is  simpler  than  the  digitizing 
system  of  Peters  and  Burke,  op.  cit,  and  includes  variable  stand 
density. 

Carter  et  al .  [23]  approached  the  optimum  road  layout 
problem  using  calculus  on  a  fuller  treatment  of  road  building  costs 
than  is  used  in  the  road  and  landing  spacing  studies  of  the  previously 
mentioned  authors. 

2.2.3  Mathematical  Programming 

In  the  last  15  years,  operations  research  tools  have  been 
applied  to  logging  problems  by  many  researchers.  Harvesting  and 
regeneration  activities  were  scheduled  using  a  linear  programming 
model  by  Curtis  [32],  but  logging  entered  the  calculations  only  as 
costs . 

Donnelly  [38]  used  linear  programming  to  determine  the 
optimum  combination  of  machines  subject  to  restraints  on  such  factors 
as  volume,  species  composition,  labour  supply  and  equipment:. 


. 

. 

4 


Lonner  T9Q]  used  linear  programming  to  minimize  storage 
and  transportation  costs  for  one  year  plans  of  wood  transportation 
activities.  A  heuristic  model  is  used  to  allocate  trucks  and  loaders 
to  individual  landings. 

More  recently,  Newnham  [118;  119]  developed  a  computer  model 
for  preparing  annual  harvesting  plans.  Linear  programming  is  used  to 
allocate  constrained  resources  of  wood  and  machines,  so  that  mill 
demands  are  met  throughout  the  year  at  minimum  cost. 

A  two-part  methodology  for  planning  cable  harvesting  areas 
and  assigning  equipment  to  them  was  developed  by  Dykstra  [40].  Stand, 
topographic  and  environmental  restraints  are  combined  with  the  mechani 
of  cable  yarding  systems  to  develop  feasible  alternatives  and  the 
harvesting  costs  for  each  block  in  a  planning  area.  An  integer 
programming  algorithm  is  used  to  design  optimum  blocks  and  to  allocate 
yarding  equipment  to  minimize  costs. 

The  main  problems  with  linear  programming  for  harvest 
planning  are  that: 

a)  variables  are  assumed  to  be  continuous, 

b)  machine  interactions  are  not  considered, 

c)  tableaux  small  enough  for  convenient  solution  but  large 
enough  for  accuracy  are  difficult  to  formulate,  and 

d)  variables  must  be  deterministic. 

Various  mathematical  programming  methods  i.e.  integer,  and  non-linear 
programming,  can  lessen  the  importance  of  the  first  three  problems, 
but  have  not  been  widely  used.  The  latter  problem  continues  to 
restrict  use  of  this  technique. 
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2.2.4  Simulation 


The  inability  of  mathematical  programming  techiniques  to 
describe  the  stochastic  relationships  common  in  harvesting  operations 
is  probably  the  greatest  single  cause  of  the  rise  simulation  as  the 
most  popular  operations  research  technique  for  investigating  harvesting 
machines  and  systems.  However,  applications  of  these  techniques  have 
also  had  problems.  In  many  studies,  the  scope  of  the  system  investigated 
was  either  too  restricted  or  too  general  to  permit  proper  evaluation  of 
the  total  harvesting  system. 

For  example,  a  simulation  model  presented  by  Gillam  [61]  models 
the  movement  of  wood  through  a  watershed,  but  does  not  discuss  how  the 
wood  gets  to  the  water.  Other  simulation  models  such  as  those  by  Leaf 
and  Alexander  [83],  and  by  Bare  [5],  treat  harvesting  as  only  one  cost 
variable  in  a  comprehensive  forest  management  model. 

The  complete  harvesting  system  model  was  not  modelled  because 
of  the  nature  of  many  harvesting  operations  and  the  overall  study 
objectives.  On  many  harvesting  operations,  the  storage  interval  at 
the  primary  landing  is  sufficiently  long  to  justify  the  separation  of 
the  harvesting  system  into  subsystems  without  invalidating  the  model. 

In  large  scale  models,  harvesting  is  only  one  activity  in  the  forest’s 
growth  cycle. 

Corcoran’s  paper  [26]  was  one  of  the  first  to  describe  a 
complete  harvesting  system  simulation  model.  However,  the  model  was 
prepared  to  illustrate  the  suitability  of  the  simulation  language 
GPSS/360,  and  does  not  appear  to  have  been  utilised  for  other  purposes. 

GASP  II  was  used  as  the  programming  language  by  Johnson  and 
Biller  [11;  70]  to  model  several  alternative  systems.  More  significant 
however,  is  that  the  model  was  validated  by  comparing  predicted  and 
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observed  results.  Like  most  simulation  models,  time  elements  were 
described  using  statistical  distributions. 

Martin’s  Timber  Harvesting  and  Transport  Simulator  differs 
in  that  regression  equations  in  conjunction  with  their  standard  errors 
were  used  to  predict  cycle  times  [96].  No  special  simulation  language 
was  used. 

The  flexibility  of  all  these  models  is  limited.  A  more 
general  model  was  developed  by  the  American  Pulpwood  Association's 
Harvesting  Research  Project  [152],  and  subsequently  expanded  by 
researchers  at  Virginia  Polytechnical  Institute  and  State  University 
[120].  Six  different  systems  can  be  simulated  by  this  model  in  five 
different  stands.  The  intent  of  the  model  was  to  evaluate  the  effects 
of  the  systems’  components  balances  and  stand  conditions  on  overall 
productivity.  Stand  conditions  were  not  intrinsic  however,  since 
production  rates  are  fixed  outside  the  model. 

Bonita  [13]  developed  a  simulation  model  for  coastal  British 
Columbia.  His  model,  designed  for  evaluating  and  comparing  existing 
or  potential  operating  policies  and  determining  equipment  requirements, 
was  validated. 

2.2.5  Network  Analysis 

Mandt  [94;  95]  determined  the  shortest  route  between  two 
points  in  a  road  network  using  a  hand  calculation  procedure  for  network 
analysis.  Wood  from  different  cutting  areas  was  allocated  between 
alternative  mills  by  minimizing  the  travel  time  between  the  cutting 


area  and  the  mill. 
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Ransing  [136;  137]  discussed  the  use  of  critical  path 
analysis  in  the  planning  of  a  forest  road  construction  project. 

Martin  [97]  designed  a  computer  program  for  analyzing  PERT  networks 
which  could  be  used  for  planning  such  projects. 

A  report  by  Carson  and  Dykstra  [22]  describes  several 
computer  programs  for  finding  the  shortest  path  through  a  transportation 
network  from  a  landing  to  a  destination.  A  programmable  desktop 
calculator  with  digitizer  was  used  to  develop  the  network.  An  algorithm 
which  finds  the  shortest  path  from  any  initial  node  to  all  other  nodes 
is  used  to  determine  the  best  transportation  route. 

Applications  of  network  analysis  have  been  limited  to  areas 
related  to  planning  and  preparing  for  the  actual  harvest  rather  than 
analyzing  alternative  systems. 

2.2.6  Production  Functions 

Multiple  regression  equations  for  various  system  components 
are  converted  to  a  probability  density  function  for  the  complete  system 
in  a  model  developed  by  Watson  and  Matthes  [155].  The  model  is  based 
on  the  assumption  that  the  system’s  productivity  equals  the  minimum 
productivity  of  the  components  for  any  level  of  production. 

2.2.7  Decision  Trees 

Woodland  [157]  used  a  decision  tree  as  a  machine  costing 
model.  He  proposed  that  probability  estimates  for  different  events 
be  prepared  and  probability  weighted  "expected"  values  be  used  in 
return  on  investment  calculations  and  risk  analysis.  The  decision  to 
purchase  a  machine  or  system  was  to  be  made  based  on  these  expected 


returns. 
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2.2.8  Others 


A  heuristic  algorithm  was  developed  by  Carlsson  [21]  to 
prepare  logging  plans  over  different  planning  horizons  in  Sweden. 

The  short-term  plans  are  basically  cutting  schedules  that  maximize 
the  utilization  of  available  machines  and  labour.  Areas  with 
constraints  on  logging  activities  are  scheduled  for  cutting  when 
the  constraint's  effects  can  be  minimized  for  the  one-year  plan. 
Individual  areas  are  allocated  to  specific  years  in  the  5-year 
plans . 

Gibson  and  Egging  [59;  41]  formulated  models  for  selecting 
the  optimal  landing  location  from  several  specified  alternatives. 
Topographic  influences  on  the  production  function  and  physiographical 
constraints  are  incorporated  in  the  model .  But  the  unit  centroids 
used  to  represent  the  cut-blocks  are  inaccurate  representatives, 
since  all  wood  from  each  block  is  assumed  to  be  transported  from 
the  block's  geometric  centre  to  the  landing.  This  appears  not  to 
recognize  that  features  which  are  constraining  at  the  centre,  may 
not  apply  to  the  complete  block,  thus  raising  average  skidding 


distance  and  costs. 
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2.3  Designing  Logging  Machines  and  Systems 

Several  approaches  have  been  applied  to  the  design  of 
machines  and  systems.  Some  of  these  approaches,  particularly 
reliability  theory,  are  exclusively  design  approaches.  However,  the 
differentiation  between  planning  and  control  models,  and  those  applied 
to  design  is  imprecise,  expecially  with  stochastic  simulation  models. 
Frequently  these  models  can  be  and  have  been  used  for  both  purposes. 
Thus  for  the  purposes  of  this  discussion,  an  attempt  has  been  made  to 
differentiate  between  each  model's  primary  objectives  while  recognizing 
the  model  may  have  other  uses. 

2.3.1  Reliability  Theory 

Two  authors  have  used  reliability  theory  to  justify  their 
approaches  to  logging  equipment  design.  Kurelek's  papers  [78;  79] 
compare  several  multi-function  machines  with  systems  composed  of 
single-function  machines.  The  improved  man-day  productivity,  lack 
of  interdependence  between  machines,  better  materials  handling 
characteristics,  and  the  redundancy  of  some  components  led  Kurelek 
to  conclude  that  multi-function  machines  have  inherent  productivity 
and  cost  advantages  over  single-function  machines. 

On  the  other  hand,  Mellgren  [105]  concluded  that  there  are 
no  significant  differences  in  the  potential  wood  production  costs 
between  different  machine  systems.  He  also  suggested  that  single¬ 
function  machines  are  more  competitive  under  adverse  terrain  conditions 
and  in  stands  of  large  trees  since  machine  size  becomes  a  limiting 


factor. 
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The  conflict  between  these  authors’  conclusions  arises 
because  of  their  treatment  of  machine  availability.  Kurelek  and 
Mellgren  generally  agree  on  the  reliability  of  different  machine 
types,  but  Kurelek' s  cost  calculations  do  not  use  these  reliability 
figures.  Both  authors  emphasize  that  the  competence  of  the  mechanical 
support  organization  controls  whether  or  not  an  operation  can  be 
successfully  mechanized. 

A  more  recent  paper  by  Mellgren  [106]  suggests  that  improved 
reliability  is  necessary  in  multi-function  machines  if  the  logging 
industry  is  to  profit  from  their  advantages.  His  results  suggest 
that  machines  with  two  or  three  functions  provide  the  best  trade-off 
between  the  advantages  of  multi-function  equipment  and  the  problem 
of  reliability  as  machine  complexity  increases. 

2.3.2  Systems  Comparison  Methodologies 

A  methodology  for  analyzing  logging  systems  was  developed 
at  Skogsarbeten  in  the  mid-sixties  [64].  The  time  required  to  complete 
various  elements  is  estimated  using  either  time  study  and/or  design 
specifications,  summed  and  converted  to  hourly  production  estimates. 

A  more  complete  machine  classification  and  expansion  of  this  methodology 
was  prepared  by  McCraw  and  Silversides  [102]  for  use  in  North  American 
conditions.  The  effects  of  operators  on  productivity  through  personnel 
delays  and  training,  and  the  interference  between  machines  are  not 
included  in  these  productivity  estimates. 


. 


- 


Boyd  and  Novak  [15]  used  the  competence  of  the  support 
organization  as  the  basis  of  a  procedure  for  evaluating  machine 
concepts  and  systems.  Eighty-four  logging  systems  were  compared  at 
four  assumed  levels  of  performance  which  reflected  organizational 
competence.  However,  this  approach  assumed  that  all  machines  are 
susceptible  to  the  same  variability  in  performance  and  availability. 
2.3.3  Simulation 

Both  deterministic  and  stochastic  simulation  models  have 
been  used  for  designing  harvesting  machines  and  systems.  Levesque  [86 
developed  a  deterministic  model  of  a  logging  truck  based  on  Newtonian 
physics  relationships  to  select  truck  components  and  to  determine 
travel  time  estimates.  Garner  and  Cooper  [58]  used  a  simular  model 
based  on  rotational  rather  than  distance  relationships  to  investigate 
the  effects  of  reducing  the  aerodynamic  drag  of  logging  trucks. 

The  Vehicle  Mission  Simulator  described  by  Gustafson  and 
Schneck  [62]  is  another  deterministic  model  used  for  design  purposes. 
Vehicle  speeds  for  specified  load,  soil  and  slope  conditions  are 
calculated  for  the  vehicle  by  solving  the  force  and  moment  equations 
[124],  The  model  is  used  to  select  optimum  components  such  as 
power-trains  and  tires  when  designing  machines. 

Fiske  and  Fridley  [44]  developed  a  model  for  selecting 
skidding  equipment  based  on  traction  characteristics.  The  model 
selects  the  skidder  or  crawler  tractor  which  produces  the  lowest 
costs  for  specific  terrain,  load  sizes  and  skidding  distances. 
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Newnham  [115;  116;  117]  has  developed  several  models  for 
use  in  designing  better  harvesting  machinery.  The  harvesting  of 
individual  trees  is  modelled  to  estimate  productivity.  However, 
the  detailed  stand  data  required  for  the  model  was  not  commonly 
available,  limiting  the  models'  potential  usefulness.  Terrain 
effects  on  machine  behaviour  were  not  modelled. 

Other  models  used  a  combination  of  stochastic  and  deter¬ 
ministic  relationships.  Damon  and  Johnson  [33]  adapted  earlier  work 
by  Almquist  [3]  in  the  modelling  of  logging  machines.  Probability 
distributions  were  used  to  generate  the  sizes  and  locations  of  trees. 
Time  elements  were  calculated  using  deterministic  relationships  on 
such  factors  as  boom  length  and  extension  speed.  The  model  was  used 
to  determine  the  optimum  specifications  for  a  thinning  machine. 

Winsauer  and  Underwood  [156]  described  a  stochastic  model 
for  estimating  the  productivity  of  a  "topwood”  harvester.  Unlike 
most  of  the  design  models,  results  of  the  simulation  runs  were 
validated  and  found  to  be  in  close  agreement  with  observed  results. 

Several  authors  have  attacked  the  problems  of  subsystem 
and  system  design  using  simulation.  Routhier's  stochastic  model  of 
the  trucking  subsystem  [140]  is  suitable  for  both  design  and  planning 
purposes.  Garner  [57]  validated  the  model  and  illustrated  its  use 
with  results  from  several  applications. 

Larsson  [80]  discussed  a  simulation  model  for  analyzing 
alternative  trucking  systems  for  design  purposes.  Regression  equations 
from  time  studies  are  used  to  estimate  elemental  times. 
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A  simulation  model  for  evaluating  a  hot  logging  system  using 
shortwood  harvesters  and  tractor-trailers  was  developed  by  Newman  [113]. 
Results  indicated  that  savings  in  delivered  wood  costs  were  possible 
over  cold  logging  systems,  but  no  validation  of  results  was  reported. 

Chip  harvesting  systems  have  been  simulated  by  Johnson  and 
Biller  [69],  and  by  Bradley  et  al .  [16;  17;  18],  Johnson’s  and  Biller's 
SAPLOS  model  was  developed  for  use  in  determining  the  best  equipment  mix 
for  a  chipping  operation.  Bradley's  simulator,  which  was  validated,  was 
designed  to  estimate  harvesting  costs  under  a  range  of  stand  conditions 
and  to  determine  the  best  equipment  mix.  Both  models  are  stochastic. 

2.3.4  Queueing  Theory 

Analytical  methods  based  on  the  idle  time  distributions  of 
machine  components  in  series  were  used  by  Meng  [107]  as  the  basis  of 
methodology  for  evaluating  design  changes.  It  is  assumed  that  both 
components  are  operating  simultaneously. 

2.3.5  Network  Analysis 

An  adaption,  of  network  analysis,  known  as  line  balancing 
was  used  by  Corcoran  [25]  to  determine  the  optimum  balance  of  components 
in  a  harvesting  system. 
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2.4  Synthesis 

The  literature  survey  indicates  that  many  analytical  techniques 
have  been  used  to  study  foresting  harvesting  machines  and  systems. 
Simulation  has  become  the  most  widely  discussed  analytical  technique. 
However,  widespread  use  had  not  followed,  suggesting  that  simulation's 
primary  disadvantages,  model  complexity  and  inability  to  determine 
optimal  solutions  directly,  are  restraining  for  potential  users. 

Most  researchers  recognized  the  interdependence  of  all 
harvesting  activities,  but  subsystems  have  generally  been  modelled. 

Such  an  approach  would  be  reasonable  if  the  model  was  intended  for  use 
in  controlling  existing  harvesting  systems.  However,  most  models  were 
developed  to  assist  in  systems  design  and  operations  planning,  tasks 
which  demand  consideration  of  many  alternatives. 

Since  none  of  the  techniques  applied  to  date  appear  to  have 
met  the  needs  of  woodlands  management,  trial  of  other  techniques  appear 
justified. 
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3.  INVESTIGATIVE  METHODOLOGY 


The  basic  problem  which  this  research  considers  is  that 
theoretically  every  forest  harvesting  system  operates  under  different 
physical,  economic  and  social  conditions.  Thus  conceivably,  there 
could  be  a  very  large  number  of  system  designs,  each  of  which  might  be 
optimal  for  some  specific  combination  of  environmental  conditions. 

Network  diagrams  which  graphically  represent  the  sequence 
of  events  iixvolved  in  completing  some  task  can  be  used  to  describe 
alternative  methods.  Various  adaptations  of  this  approach  have  been 
used  to  design  agricultural  systems,  but  no  reference  was  found  to  its 
application  to  selecting  forest  harvesting  systems. 

Preston  [134;  135]  adapted  the  critical  path  method  (CPM) 
so  that  both  time  and  cost  variables  were  evaluated  simultaneously 
in  applications  to  the  layout  and  operation  of  piggeries  and  to  select 
optimal  irrigation  systems.  Coupland  and  Halyk  [30]  used  network 
analysis  to  evaluate  forage  harvesting  systems.  A  simulation  model 
was  used  to  select  alternative  paths  through  a  network  of  forage 
handling  alternatives  by  Lievers  [87], 

A  computer  program  (SPNA)  initially  developed  by  Preston  [133] 
for  selecting  the  shortest  path  through  a  network,  was  adapted  by  Lievers 
[88]  to  rank  all  paths  through  the  network.  Ogilvie  et  al .  [121]  used 
SPNA  to  evaluate  manure  handling  systems  for  swine  and  dairy  cattle. 

Safley  and  Price  [141],  and  Burney  et  al .  [20]  have  used 
different  versions  of  the  network  approach  to  analyze  alternative 


manure  handling  systems. 
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The  shortest  path  network  approach  to  system  design  was 
selected  for  this  research  because  of  the  following  advantages  apparent 
from  the  literature: 

a)  alternative  combinations  and  sequences  of  equipment 
can  be  easily  outlined  in  network  form, 

b)  many  alternatives  can  be  quickly  and  cheaply  evaluated, 

c)  optimal  and  near-optimal  system  designs  are  identified 
for  more  intensive  study  if  desired, 

d)  non-linear  functions  can  be  used,  and 

e)  sensitivity  analysis  can  be  readily  undertaken. 
Alternative  approaches  such  as  simulation  and  mathematical  programming 
do  not  have  all  these  advantages. 

The  methodology  which  follows  uses  shortest  path  network 
analysis  to  assist  the  forest  harvesting  system  designer  identify  the 
optimal  system  for  his  conditions.  Alternative  combinations  of 
logging  equipment  are  compared  on  a  common  basis  to  permit  the 
selection  of  the  least-cost  system  or  systems.  The  procedures  used 
will  be  described  under  the  following  subsections: 

3.1  Preparation  of  the  Network 

3.2  Establishment  of  Activity  Durations 

3.3  Shortest  Paths  Network  Analysis  Program 
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3.1  Preparation  of  the  Network 

A  network  is  a  graphical  model  which  defines  the  sequence 
of  events  and  activities  between  some  initial  state  and  some  concluding 
state.  An  event  is  the  identifiable  point  at  which  some  defined  state 
exists,  and  is  denoted  on  a  network  diagram  by  a  circle  or  "sausage", 
called  a  node.  An  activity,  represented  on  the  graph  by  an  arrow,  is 
an  operation  which  is  required  to  reach  some  defined  state,  and  is 
called  a  branch,  an  arc  or  a  link.  An  activity  having  no  quantity 
associated  with  it  but  needed  to  maintain  a  logical  sequence,  is 
represented  by  a  dashed  arrow  called  a  dummy.  Dummies  are  also 
required  for  uniqueness  and  to  avoid  ambiguity  when  there  are  several 
alternative  activities  between  two  nodes.  The  network  is  the  graph 
produced  by  connecting  sequential  events  with  arcs  representing  the 
activities  involved  in  producing  the  change  in  state. 

Figure  6  illustrates  a  generalized  network  diagram  for 
alternative  tree-length  harvesting  systems.  The  initial  state  in 
the  tree-length  system  was  assumed  to  be  the  tree  standing  in  the 
forest.  Several  alternative  operations  and  transportation  activities 
convert  from  a  standing  tree  into  a  tree-length  log,  transferring  it 
from  the  stump  to  some  intermediate  landing  and  thence  onto  the 
processing  plant.  Along  the  way,  it  may  or  may  not  be  temporarily 
stored.  The  unloaded  tree  at  the  mill  represents  the  final  state 


and  node. 
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Figure  6.  General  network  diagram  of  alternative  activities  in  tree-length  harvesting  systems. 
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Because  of  network  size  and  the  mass  of  data  required  to 
describe  all  alternative  systems,  this  research  was  limited  to  a 
tree-length  harvesting  system  as  it  could  be  applied  in  Eastern  Canada. 
As  such,  the  "optimal"  system  in  the  remainder  of  this  thesis  is  the 
least-cost  tree-length  harvesting  system  which  might  not  be  optimal 
if  all  alternative  systems  were  considered.  However,  shortest  path 
network  analysis  is  a  general  technique  and  can  be  easily  adapted  for 
use  with  different  networks  and  inputs  for  other  systems. 

Figure  7  shows  an  expansion  of  the  general  network  diagram 
outlining  the  alternative  machines  and  subsystems  which  can  perform 
the.  various  operation  and  transport  activities  of  the  tree-length 
system. 

Three  assumptions  were  made  to  reduce  the  network's  complexity 
while  preparing  the  expanded  network: 

a)  truck  transportation  is  the  only  method  available  for 
moving  the  wood  from  landing  to  mill.  Since  the  mode 
of  transportation  is  usually  predetermined  by  the 
availability  of  alternatives  e . g ,  railroads,  rivers, 
and  road  systems,  and  most  wood  is  trucked  for  some 
distance  in  Canada,  this  is  a  reasonable  simplification; 

b)  hot-logging  systems,  where  some  equipment  would  off-load 
onto  trucks  are  infeasible.  Interference  problems  between 
machines  have  prohibited  use  of  hot-logging  systems  in 


most  areas;  and 
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Figure  7.  Expanded  network  diagram  of  alternative  tree-length  harvesting  systems. 
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c)  temporary  storage  activities  usually  have  durations  less 
than  a  week  so  that  no  interest  charges  are  added  to 


42 


production  costs  in  lieu  of  the  capital  tied-up  in 
uncompleted  products.  Lengthy  storage  between  activities 
is  common  only  between  the  harvesting  and  transportation 
subsystems.  This  delay  varies  from  a  few  days  to  as 
much  as  a  year  depending  on  the  specific  operation. 

All  systems  are  likely  to  have  a  delay  of  similar  duration 


at  this  point. 
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3.2  Establishment  of  Activity  Durations 

Each  arc  in  the  network  diagram  of  the  alternative  systems 
shown  in  Figure  7  was  assigned  a  duration.  In  this  analysis,  activity 
durations  were  measured  as  wood  cost  (CW)  expressed  as  a  cost  per  cubic 
metre  of  wood  produced. 

The  general  relationship  for  the  duration  of  any  activity 
reflects  the  economics  and  productivity  of  the  machinery  and  labour 
involved  in  the  activity. 


_  HCOST 
HP  ROD 


(1) 


where:  CW  =  wood  cost  for  the  activity  ($/m3) 

HCOST  =  hourly  costs  for  the  machinery  and  labour  used  in 
the  activity  ($/productive  machine  hour) 

HPROD  =  hourly  productivity  of  the  machinery  and  labour  used 
in  the  activity  (m3 /productive  machine  hour) . 

The  derivation  of  the  equations  used  to  estimate  the  hourly  costs  and 
productivity  for  the  activities  are  discussed  in  the  next  two  sections. 
The  duration  of  each  activity  measured  in  dollars  per  cubic  metre  of 
wood  produced  were  calculated  using  equation  1  from  the  estimates  of 
hourly  costs  and  productivity. 
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3.2.1  Cost  Prediction  Equations 

A  costing  method  based  on  optimal  economic  machine  life  was 
used  to  predict  equipment  costs  on  a  before-tax  basis.  This  approach 
was  selected  because  the  tendency  for  maintenance  costs  to  increase  and 
for  productivity  to  decrease  as  a  machine  ages,  results  in  an  economic 
life  that  is  often  shorter  than  the  physical  life. 

The  prediction  of  logging  equipment  costs  considered  the 
following  items: 

a)  purchase  price, 

b)  machine  life, 

c)  salvage  value, 

d)  cost  of  capital  (interest  rate) , 

e)  insurance  costs, 

f)  licensing  costs, 

g)  labour  costs, 

h)  fuel,  lubricant  and  hydraulic  oil  costs, 

i)  maintenance  costs,  and 

j)  machine  utilization. 

Costing  models  used  previously  in  studies  of  logging  machines 
and  systems  have  faults  related  to  their  treatment  of  capital  costs. 

The  straight  write-off  approach  used  by  Boyd  and  Novak  [15]  makes  no 
charge  for  invested  capital  and  assumes  no  salvage  value.  The.  method 
used  in  PAPRICAN’s  and  FERIC’s  machine  evaluation  studies  e .g.  [128;  85], 
write-off  plus  interest  on  half  the  investment,  assumes  no  salvage  value 
and  equal  annual  interest  charges.  The  widely  used  straight-line  write¬ 
off  plus  interest  on  the  average  investment  e. g.  [6;  24;  154]  fails  to 
recognize  compound  interest  effects. 


All  these  methods  assume  an  equal 
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annual  write-off  fixed  by  an  estimate  of  recovery  period  which  is  usually 
defined  as  the  machine's  "useful"  life. 

This  research  used  a  capital  recovery  with  return  approach 
to  calculate  capital  costs.  In  comparison,  the  straight-line  write-off 
plus  interest,  on  the  average  investment,  the  most  accurate  of  the 
previously  discussed  costing  methods,  provides  capital  cost  estimates 
which  are  about  11%  low  for  a  10  year  recovery  period  and  a  12%  interest 
rate.  The  capital  costs  (CC)  expressed  In  annual  equivalent  dollars 
per  productive  machine  hour  (PMH)  were  calculated  using  the  formula 


(PP-SV.)  CRF  +  SV.  x  I 

3 


SHY  x  U 


(2) 


\tfhere :  CC . 

3 

PP 

SV. 

3 

CRF 


capital  costs  in  year  "j " 
purchase  price  ($) 
salvage  value  in  year  "j " 
capital-recovery  factor  = 


($/PMH) 


1(1  +  I)3 

(1  +  I)3-l 


I  =  interest  rate  (decimal  %) 

SHY  =  scheduled  hours  yearly  (scheduled  machine  hours) 

U  =  utilization  (decimal  %) . 

A  machine's  salvage  value  is  a  measure  of  the  remaining 
utility  in  the  machine.  Value  is  lost  over  time  as  the  efficiency  of 
the  machine  decreases.  For  example,  an  engine  reaches  its  maximum 
efficiency  following  a  short  run-in  period  after  which  its  overall 
efficiency  declines  as  the  cylinder  walls,  pistons  and  rings  wear. 

The  machine  always  maintains  some  value  if  only  as  scrap  metal  though 
the  costs  of  disposal  may  negate  any  salvage  value. 


, 


For  this  analysis,  a  machine's  value  was  assumed  to  decrease 
at  an  annual  rate  of  30%.  The  salvage  value  of  a  machine  in  year  " j " 
(SV  )  was  es  timated  to  be 


SV .  =  PP  (1  -  0.30)- 


(3) 


The  remaining  costs  of  harvesting  equipment  consists  of  both 
fixed  and  variable  cost  items.  Licensing,  insurance  and  to  a  lesser 
extent  operator  wages  are  fixed  and  are  incurred  whether  or  not  the 
machine  operates.  The  variable  costs  of  fuel,  oils  and  lubricants, 
and  repairs  are  dependent  on  the  number  of  hours  worked  by  the  machine 
All  these  costs  were  converted  to  a  productive  machine  hour  basis. 

Insurance  costs  (INC) ,  an  expenditure  for  protection  against 
damage  to  the  machine,  were  assumed  to  be  a  constant  3%  of  the  machine 
residual  value  annually. 


INC 


0.03  x  SV. 
_ 1 

SHY  x  U 


(4) 


Licencing  costs  (LC)  were  included  only  for  the  truck  and 
trailer  in  this  analysis.  This  cost  was  fixed  at  $1700  per  annum 
based  on  the  vehicle  registration  fees  of  the  province  of  Ontario. 
Other  equipment  would  not  usually  be  licensed  for  highway  travel. 


LC 


$1700 


(5) 


SHY  x  U 
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Both  operator  wages  and  fringe  benefits  based  on  a  recent 
union  contract  [4]  were  included  in  manpower  costs  (MC) .  A  fringe 
benefit  factor  (FBF)  equalling  35%  of  direct  wages  was  applied.  It 
was  assumed  that  manpower  costs  would,  continue  to  rise  at  an  average 
rate  of  5%  annually  in  constant  dollars  as  they  have  in  the  past 
twenty  years.  Hourly  manpower  costs  ($/PMH)  were  calculated  using 
the  equation 

„„  _  W  (1  +  FBF)  x  a.05)3'1 

_  -  (b) 

where:  W  =  operator's  wage  ($/scheduled  machine  hour). 

The  cost  of  fuel  (FC)  is  dependent  on  engine  design 
characteristics ,  speed  and  efficiency,  fuel  density,  fuel  price  and 
the  load  factor  (LF) .  The  load  factor,  a  measure  of  the  severity  of 
the  job  and  conditions  under  which  the  machine  is  operating,  is  defined 
as  the  actual  fuel  consumption  as  a  proportion  of  the  maximum  capacity 
of  the  engine  to  burn  fuel.  Sample  load  factors  typical  of  harvesting 
equipment  are  summarized  in  Table  3.  Hourly  fuel  costs  ($/PMH)  were 
estimated  using  the  following  equation. 

FC  =  0.21  x  GP  x  LF  x  FP  (7) 

where:  GP  =  gross  engine  power  (kW) 

FP  =  fuel  price  ($/litre). 

It  was  assumed  that  fuel  consumption  occurs  only  when  the 
machine  is  producing  wood.  In  reality,  some  non-productive  operating 
time  does  occur,  bx.it  it  usually  represents  less  than  3%  of  scheduled 
time.  As  such,  its  effect  on  costs  is  insignificant  and  affects  all. 


machines  approximately  equally. 
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The  cost  of  oil.  and  lubricants  (OLC)  varies  with  engine 
power  and  the  capacity  and  complexity  of  the  hydraulic  system.  These 
costs  are  commonly  expressed  as  a  proportion  of  fuel  costs  since  a 
strong  relationship  exists  between  the  unit  prices  of  fuel,  oils  and 
lubricants . 

OLC  =  FC  x  HCF  (8) 

where:  HCF  =  hydraulic  complexity  factor  (see  Table  3). 

The  average  hourly  repair  costs  (RC)  over  a  machine’s  life 
were  estimated  using  the  following  equation  [24,  p.  28-29;  15,  p.  5], 


RC 


PF  x  RCF 
1000 


(9) 


where:  RCF  =  repair  cost  factor  ($/$1000  of  purchase  price/PMH) . 

This  cost  includes  all  parts  and  mechanics’  labour.  The  repair  cost 
factor  depends  on  machine  type,  operating  conditions  and  the  competence 
of  the  maintenance  organization.  Typical  values  for  this  factor  are 
shown  in  Table  3. 

A  constant  repair  cost  factor  was  used  because  long-term 
patterns  of  repair  costs  for  logging  equipment  tend  to  be  constant 
after  adjustment  for  inflation  until  the  frame  fails.  Boyd  suggests 
that  the  reduced  repair  time  with  age  results  from  corrective  action 
on  problem  areas  and  the  build-up  of  maintenance  experience  [14,  p.  17]. 
In  addition,  much  of  the  machine  is  not  as  old  as  the  chassis  since 
components  are  continually  replaced  during  repairs. 

Total  hourly  operating  costs  (OC)  are  calculated  by  summing 
insurance,  licensing,  manpower,  fuel,  oil  and  lubricant  and  repair  costs. 


OC  =  INC  +  LC  +  MC  +  FC  +  OLC  +  RC 


(10) 
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Utilization,  that  is  the  proportion  of  scheduled  hours  that 
the  machine  is  actually  productive,  is  dependent  upon  the  mechanical 
reliability  inherent  within  the  machine,  and  the  operational  effectiveness 
of  the  woodlands  organization. 

One  measure  of  a  machine rs  mechanical  reliability  is  its 
availability,  that  is  the  portion  of  total  scheduled  time  in  which  the 
machine's  mechanical  condition  permits  work.  The  availability  (A)  of 
a  machine  can  be  defined  as 


MTBR 

MTBR  +  MTTR 


(11) 


where:  MTBR  =  mean  time  between  repairs 
MTTR  =  mean  time  to  repair. 

Values  used  in  this  research  for  the  MTBR  of  logging  equipment,  were 
gathered  from  various  long-term  studies  of  machine  repair  statistics 
e. g.  [51;  65;  132].  Where  values  were  unavailable,  estimates  were 
made  by  combining  available  information  for  similar  components  on 
other  equipment  [78;  106]. 

A  mean  time  to  repair  (MTTR)  of  1.8  hours  was  assumed  in 
determining  the  mechanical  availability  [14,  p.  17].  The  MTTR  varies 
with  machine  and  mechanical  characteristics  such  as  accessibility  and 
part  complexity,  the  mechanic’s  skill  and  the  facilities  with  which  he 
works,  but  the  1.8  hour  value  is  typical  of  most  harvesting  machines 
and  woodlands  maintenance  organizations . 

Availability  is  not  directly  convertible  to  utilization, 
since  some  repairs  may  be  completed  outside  scheduled  hours.  Some 
companies’  maintenance  policies  restrict  repairs  to  scheduled  hours, 
but  in  many  cases  85%  of  all  repair  time  occurs  within  scheduled 
hours  [48,  p.  4;  49,  p.  4]. 


.  ■ 


.r  ' 


51 


In  addition  servicing  activities  such  as  fuelling,  cleaning, 
oil  changes  and  preventive  maintenance,  are  required  on  the  machine. 

The  time  required  for  these  activities  varies  with  the  manufacturer's 
recommended  servicing  schedule  and  on  company  policy,  but  is  proportional 
to  the  machine's  productive  hours.  Typically,  10%  of  productive  hours  are 
spent  servicing  the  machine.  Availability  thus  becomes 


_ MTBR _ 

1.1  x  MTBR  +  0.85  x  MTTR 


(12) 


Operational  effectiveness  (OE)  is  the  other  factor  which 
affects  the  machine  utilization.  This  factor  measures  the  organizational 
and  human  controlled  performance  of  the  machine. 


OE  = 


PMH 


PMH  +  WM  +  WP  +  NMD 


(13) 


where:  WM  =  time  spent  waiting  for  mechanic (s)  or  repair  facilities, 
WP  -  time  spent  waiting  for  repair  parts, 

NMD  =  non-mechanical  time  delay. 

Boyd  [14,  pp.  50-51]  suggested  that  non-mechanical  delay  time 
could  be  divided  into: 

a)  Operational  lost  time,  caused  by  weather  or  terrain 
conditions  (e. g .  warm-up  or  machine  stuck) ,  interference 
between  other  machines  in  the  system,  waiting  for 
supervisor's  instructions,  aiding  other  machines,  etc . ; 

b)  Personal  (e. g .  operator  late,  sick,  rest  allowances);  and 

c.)  In-shift  moving  (e.  g.  moving  between  cut  blocks)  . 


. 

. 
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Various  studies  of  operational  effectiveness  have  shown  it  to  be 
extremely  variable  ranging  from  54%  [14,  p.  26]  to  95%  [65,  p.  4]. 
Values  of  operational  efficiency  appropriate  to  each  machine  were 
used  in  the  analysis. 

The  product  of  availability  and  operational  effectiveness 
is  utilization  (U) . 


U  =  A  x  OE 


MTBR 


1.1  x  MTBR  +  0.85  x  MTTR 


x  OE 


(14) 


The  utilization  factor  was  used  to  convert  costs  from  a  scheduled  hour 
to  a  productive  hour  basis  where  necessary. 

The  sum  of  the  capital  and  operating  costs  are  the  total 
machine  costs.  However,  since  both  vary  with  machine  age,  the  cash 
flows  were  converted  to  a  common  time  basis.  The  operating  costs  for 
each  year  are  discounted  and  summed,  and  the  total  converted  to  annual 
equivalent  dollars  per  productive  machine  hour. 


J  _v 

EAC.  =  CC.  +  CRF  x  (  l  OC.  x  (1+i)  ) 

33  k=l 

where:  EAC,  =  equivalent  annual  cost  in  year  M j M  ($/PMH) 

1 

0C1  =  operating  costs  in  year  "k"  ($/PMH) 


(15) 


Hourly  operating  costs  were  measured,  in  constant  dollars,  so  an 
"inflation-free"  interest  rate  of  12%  was  used  in  the  analysis. 


. 
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The  machine  age  which  produced  the  minimum  equivalent  periodic 
cost  for  equation  15  is  the  optimum  economic  life  [43,  pp .  125-147; 

139,  pp.  90-214].  It  was  assumed  that  obsolescence  due  to  technological 
changes  necessitate  machine  replacement  within  a  maximum  of  10  years. 

If  the  lowest  cost  point  had  not  been  reached  by  that  point,  the  optimum 
life  was  assumed  to  be  10  years.  The  machine  hourly  cost  (HCOST)  used 
to  calculate  activity  durations  (equation  1)  was  the  equivalent  periodic 
cost  for  the  machine  kept  for  its  optimal  life. 

HCOST  =  Minimum  .  _  (EAC . )  (16) 

3  -1-hLO  2 


' 

. 
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3.2.2  Productivity  Prediction  Equations 

The  productivity  of  forest  harvesting  equipment  is  dependent 
on  a  basic  cycle  time  inherent  to  the  machine  which  is  modified  by 
conditions  in  the  physical  environment  and  by  operator  performance. 

As  such,  a  model  of  productivity  is 


HPROD  = 


Q  x  6000 


where:  HPROD 

Q 

CT 


CT  (1  +  EF  +  OP) 
hourly  productivity  (m3/PMH) 
quantity  produced  per  cycle  (m:i) 
machine  cycle  time  (centiminutes  per  cycle) 


(17) 


EF  =  productivity  adjustment  factor  for  physical 
environmental  conditions 

OP  =  productivity  adjustment  factor  for  operator  performance. 
However,  the  potential  usefulness  of  this  model  is  limited. 

Machine  design  itself  fixes  the  potential  zone  of  application 
for  a  piece  of  harvesting  equipment.  For  example,  a  f eller-buncher  can 
operate  only  where  the  slope  is  less  than  some  critical  angle  above  which 
the  machine  will  roll  over.  Similarly  a  skidder  can  pull  a  larger  load 
only  while  its  tractive  effort  exceeds  the  forces  resisting  skidding. 
Since  performance  decreases  rapidly  as  these  physical  limits  are  reached, 
machines  tend  to  be  used  in  areas  where  performance  is  least  affected  by 
physical  conditions.  Within  each  zone  of  application,  tree  parameters 
such  as  stem  volume  are  frequently  the  only  influencing  factor. 

As  previously  discussed,  quantifying  environmental  influences 
has  proved  to  be  extremely  difficult  because  of  interactions  between 
various  conditions,  and  the  variability  of  the  micro-environment. 
Equipment  such  as  trucks  and  loaders  which  operate  under  relatively 


constant  conditions  are  the  only  consistent  exception. 


' 
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Quantifying  factors  affecting  operator  performance  has  proved 
equally  difficult.  While  many  authors  e . g ,  [85]  have  observed  operator 
effects  on  cycle  times,  these  differences  could  frequently  be  attributed 
to  inexperienced  operators.  The  study  of  tractor-mounted  shears  by 
Cot tell  et  al.  [29]  was  the  only  study  found  during  the  literature 
review  which  quantified  the  impact  of  factors  such  as  experience, 
motivation  and  manual  dexterity  on  operator  performance. 

The  best  available  predictive  equations  for  machine  cycle 
time  have  coefficients  of  multiple  determination  (R^)  which  explain 
less  than  65%  of  cycle  time  variation.  Frequently  only  a  small  portion 
of  the  cycle  time  was  variable,  with  machine  design  and  unmeasured  or 
unused  variables  producing  large  constants  in  the  equations.  While 
recognizing  the  problems  involved,  this  research  used  the  etvailable 
information  on  cycle  times  taken  in  most  cases  from  various  machine 
evaluation  studies. 

As  such,  a  simpler  model  of  productivity  was  used 


HP  ROD 


MV  x  6000 
CT 


(18) 


where:  MV  =  merchantable  volume  per  tree  (m3)  . 

The  equations  for  machine  cycle  time  (CT)  used  in  this  analysis  are 
summarized  in  Appendix  B.  The  hourly  productivity  (HPR.0D)  of  each 
activity  was  calculated  using  equation  18. 

Since  the  productivity  and  most  cycle  time  calculations 
required  various  tree  parameters,  relationships  for  stump  diameter, 
merchantable  tree-length,  green  weights  of  tree-lenths  and  full-trees, 
and  merchantable  volume  were  identified  based  on  tree  diameter  breast 
height  (DBH)  .  These  equations  are  summarized  in  Appendix  B.2.  for  the 
four  main  softwood  species  in  Eastern  Canada. 


. 
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3.2.3  Example  of  Activity  Duration  Calculations 

To  illustrate  how  the  hourly  cost  and  productivity  estimates 
were  related,  the  activity  duration  calculation  for  the  arc  joining 
nodes  1  and  142  in  Figure  7  follows. 

The  machine  used  in  this  activity  is  a  Koehring  KFF  feller 
forwarder.  Table  4  summarizes  the  costing  data  used  in  the  sample 
calculation.  The  general  variables  apply  to  all  activities.  Activity 
variables  are  specific  to  the  machine  used  in  the  activity. 

By  substituting  the  costing  variables  in  Table  4  into  equations 
2-14,  the  hourly  operating  and  capital  costs  shown  in  Table  5  were 
calculated.  Since  the  operating  costs  are  a  future  value,  they  were 
converted  to  an  annual  equivalent  before  summation  with  capital  costs 
to  obtain  total  hourly  costs  (equation  15) . 


Table  4.  Costing  Variables  for  the  Sample  Activity  Duration  Calculation. 
General  Variables : 

Interest  Rate  (%)  .  12 

Fuel  Price  ($/litre)  .  0.22 

Repair  Time  In-Shift  (%  of  total  repair  time)  .....  85 

Fringe  Benefit  Factor  (%  of  hourly  wage)  . .  35 

Insurance  Rate  Factor  (%  of  unrecovered  value)  ....  3 
Rate  of  Lost  Utility  (%  per  year)  .................  30 


Activity  Variables: 


Scheduled  Hours  Yearly  (SMB)  .  . . . .  3840 

Purchase  Price  ($)  . .  325  000 

Hourly  Wage  ($/SMH)  . 9.66 

Gross  Power  (kW)  . 257 

Load  Factor  .  0.70 

Hydraulic  Complexity  Factor . . . . .  0.25 

Repair  Cost  Factor  . . . . .  0.10 

Mean  Time  Between  Repairs  (PMH)  . . . .  11.8 

Operational  Effectiveness  . 0.90 

License  Cost  ($/yr)  . 0 
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Table  5.  Equivalent  Annual  Hourly  Cost  ($/PMH)  for  the  Sample  Calculation. 
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In  the  example,  a  minimum  cost  point  was  not  reached  within 


a  10  year  period.  Thus,  the  hourly  cost  (HCOST)  of  the  activity  was 
equated  to  the  equivalent  hourly  cost  in  year  10,  $85.61  per  productive 
machine  hour. 

Table  6  summarizes  the  operating  variables  used  to  calculate 
hourly  productivity.  The  values  for  the  weight  of  full  trees  and 
merchantable  volume  per  tree  were  calculated  for  average  tree  size 
using  the  tree  parameter  equations  in  Appendix  B.2  for  black  spruce. 


Table  6.  Productivity  Variables  for  the  Sample  Activity  Duration 


Calculation. 

Operating  Variables: 

Species  . ......... . . . . .  Black  Spruce 

Average  Tree  Size  (DBH  in  cm)  . . .  18.00 

Average  Primary  Transport  Distance  (meters)  .  2.00 

Calculated  Variables : 

Weight  of  Full  Tree  (kg)  . . .  195.6 

Merchantable  Volume  per  Tree  (m3)  .  0.152 
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The  cycle  time  per  tree  (CT)  for  the  Koehring  KFF  was 
calculated  using  the  equation  [85,  p.  5] 


CT  =  49.4  + 


(52  +  5.84  x  DIST) 

NT 


(19) 


-  7-2  x  (INT  « 


=  42.4  x  MV 


=  52.76  cmin/tree 


where:  DIST  =  one-way  travel  distance  (200  x  1.4  =  280  m) 

29000 

NT  =  number  of  trees  per  trip  =  INT  (■■----  )  =  148 

w  1  r  i 

INT  =  integer  function  which  returns  the  non-decimal 
portion  of  a  number 
WTFT  =  weight  of  full  trees  (kg) 

MV  =  merchantable  volume  per  tree  (m3)  . 

Substituting  the  cycle  time  into  equation  18  gives  the  hourly  productivity 


HP  ROD  = 


MV  x  6000 
CT 


=  0.152  x  6000 
52.76 

=  17.3  m3/‘PMH 


Combining  the  hourly  cost  and  productivity  estimates  in 
equation  1  gives  the  activity's  duration. 


CW  = 


HCOST 

HPROD 


85.61 


17.3 

=  $4. 96/m3 


; 


3.3  Shortest  Path  Network  Analysis  Program 

An  interactive  package  of  BASIC  computer  programs  entitled 
LOGNAP,  an  acronym  for  LOGing  Network  Analysis  Programs,  was  developed 
to  analyze  a  network  of  alternative  systems  such  as  that  developed  in 
Section  3.1.  The  programs  define  the  network,  calculate  the  activity 
durations  and  rank  a  user-specified  number  of  least-cost  paths  through 
the  network. 

The  flowchart  shown  in  Figure  8  summarizes  the  basic  logic 
sequence  of  LOGNAP.  The  complete  program  listing  and  some  explanation 
of  the  non-standard  features  in  the  particular  BASIC  interpreter  used 
in  the  research  are  contained  in  Appendix  C. 

The  sections  which  follow  discuss  the  limitations  on  the  size 
of  network  to  be  analyzed,  and  algorithms  used  to  renumber  nodes  and 
find  the  shortest  paths. 

3.3.1  Program  Limitations 

Five  basic  parameters  are  entered  to  describe  the  network 
and  control  processing: 

a)  the  number  of  activities  in  the  network, 

b)  the  number  of  nodes  in  the  network, 

c)  the  number  of  least-cost  paths  desired  by  the  user, 

d)  an  estimate  of  the  maximum  number  of  nodes  on  any 
least-cost  path,  and 

e)  the  number  of  sets  of  variables. 

As  presently  dimensioned,  the  number  of  activities  and  nodes  in  the 
network  is  limited  to  250,  a  maximum  of  50  paths  can  be  selected,  and 
a  path  cannot  contain  more  than  25  nodes. 
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Figure  8,  Flowchart  outlining  the  logic  sequence  of  LOGNAP 


3.3.2  Node  Renumbering  Routine 


The  subroutine  RNUM  renumbers  the  network  nodes  to  satisfy 
the  so-called  I-J  rule4.  This  was  included  in  LOGNAP  to: 

a)  simplify  the  network  development  process  by  permitting 
haphazardly  numbered  nodes  to  be  used, 

b)  facilitate  modification  of  the  network  as  required,  and 

c)  satisfy  the  requirements  of  the  shortest  path  algorithm. 
The  routine  was  one  adapted  by  Burney  [20]  from  one  developed  earlier 
by  Lievers  [88]. 

The  nodes  of  the  network  are  numbered  between  1  for  the 
origin,  that  is  the  node  having  no  predecessors,  and  the  total  number 
of  nodes  in  the  network  for  the  terminus  using  all  numbers  between 
these  values. 

The  following  are  the  main  steps  in  the  routine: 

STEP  (a)  -  identify  the  terminal  node  number  of  the  network,  and  set 
counters. 

STEP  (b)  -  increment  activity  counter,  and  test  if  all  activities  have 
been  checked.  If  so,  transfer  to  Step  (g) . 

STEP  (c)  -  identify  an  activity  having  the  terminal  node,  and  equate 
the  terminal  node  of  the  activity  to  the  highest  remaining 
terminal  node  number.  Note  that  during  the  initial  pass, 
the  highest  remaining  terminal  node  number  is  the  number 
of  nodes  in  the  network. 


The  I-J  rule  requires  that  the  nodes  in  a  network  be  numbered  so 
that  the  initial  node  number  of  any  activity  is  less  than  the 
terminal  node  of  that  activity. 
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STEP  (d) 


STEP  (e) 


STEP  (f) 


STEP  (g) 


STEP  (h) 


test  if  the  other  activities  have  the  same  initial  node. 

A  unique  initial  node  implies  that  no  alternative  path 
exists  between  the  initial  and  terminal  nodes  of  the 
activity.  If  the  node  is  not  unique,  transfer  to  Step  (b) . 
decrement  the  highest  remaining  initial  node  number,  and 
set  the  initial  nodes  of  all  activities  having  that  initial 
node  to  the  highest  remaining  node  number, 
if  all  initial  nodes  have  been  renumbered,  go  to  Step  (h) . 
If  not,  transfer  to  Step  (b)  . 

decrement  the  highest  remaining  terminal  node  number,  and 
identify  the  initial  node  of  the  highest  remaining 
uprenumbered  activity.  Reset  activity  counter  and  transfer 
to  Step  (b) . 

reset  the  initial  nodes  to  their  original  values,  and 
establish  a  key  for  the  renumbered  nodes. 
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3.3  Yen's  Algorithm 

The  shortest  paths  through  the  network  are  determined  using 
an  algorithm  developed  by  Yen  [159].  This  algorithm  increases  its 
computational  needs  linearly  with  relation  to  the  number  of  paths 
sought,  whereas  other  approaches  to  the  shortest  path  problem,  e . g . 

Lievers  [88]  increase  their  computational  requirements  exponentially. 

Yen's  algorithm  capitalizes  on  the  fact  that  once  the  shortest 
path  has  been  found,  the  nodes  on  any  sub-optimal  path  must  coincide 
with  at  least  one  node  on  the  previously  identified  paths.  Thus  to 
find  the  next  shortest  path,  only  the  shortest  deviations  from  the 
previously  found  paths  need  to  be  evaluated.  The  shortest  of  these 
candidates  is  the  next  best  path. 

As  proposed  by  Yen,  two  lists  (matrices)  are  used  to  store 
the  best  paths  found  so  far  and  the  candidate  paths  respectively.  In 
LOGNAP  a  modification  which  reduces  the  storage  requirements  has  been 
made  to  the  algorithm.  Since  only  the  desired  number  of  paths  need 
to  be  stored,  both  lists  can  be  saved  in  the  same  matrix.  The  best 
paths  are  stored  in  the  top  portion  of  the  matrix  and  the  candidates 
in  the  bottom  part.  As  each  path  is  identified,  it  is  inserted  in  the 
matrix  at  the  appropriate  position  from  lowest  to  highest  cost  and  any 
longer  paths  are  forced  down  in  the  stack  or  even  dropped.  Upon  completion 
of  each  iteration,  the  shortest  candidate  path  is  added  to  the  bottom  of 
the  permanent  list  by  incrementing  a  row  counter  by  one. 
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The  basic  steps  in  the  algorithm  are  as  follows : 

STEP  (a)  -  find  the  shortest  path  through  the  network.  In  LOGNAP 
the  shortest  path  between  any  two  nodes  was  found  using 
Elmaghraby's  algorithm  (see  next  section). 

STEP  (b)  --  for  each  subsequent  path  desired,  check  if  the  subpath 

between  each  node  of  the  most  recently  found  path  and  the 
origin  node  coincide  with  a  portion  of  any  previously  found 
paths.  If  so,  set  the  duration  of  the  activity  following 
the  match  to  a  large  number. 

STEP  (c)  -  find  the  shortest  path  around  the  blockage  to  the  terminal 
node.  This  subpath  becomes  the  spur. 

STEP  (d)  -  form  a  new  path  by  joining  the  portion  of  the  path  behind 
the  blockage,  the  root,  with  the  spur  and  add  the  path  to 
the  list  of  candidates. 

STEP  (e)  -  select  the  shortest  path  from  the  list  of  candidates  and  add 
to  the  permanent  list. 

STEP  (f)  -  continue  at  Step  (b)  until  the  number  of  desired  paths  has 
been  found. 

Satisfaction  of  the  I-J  rule  is  necessary  for  Yen's  algorithm. 
Before  the  paths  are  outputted.,  the  original  nodes  are  restored  using 
the  key  developed  in  the  node  renumbering  routine. 

The  formulation  of  Yen’s  algorithm  used  in  LOGNAP  is  based 
on  a  program  GSPNA  developed  by  Horne  [67].  Although  GSPNA  is  less 
efficient  than  a  second  program  GMSP  also  prepared  by  Horne,  storage 
limitations  made  it  preferable  on  the  computer  used  in  this  research. 

The  more  efficient  GMSP  version  of  the  algorithm  can  be  easily  adapted 
for  use  in  LOGNAP  if  run  on  a  larger  computer  (>64K  of  core  memory) . 
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3.3.4  Elmaghraby 1 s  Algorithm 

Numerous  algorithms  have  been  developed  for  finding  the 
shortest  path  between  two  specified  nodes,  e . g .  Bellman  and  Kalaba  [7], 
Dantzig  [34],  Dijkstra  [36],  Ford  and  Fulkerson  [54],  and  Minty  [109]. 

A  1969  review  of  shortest  path  algorithms  made  by  Dreyfus  [39]  concluded 
that  the  Dijkstra  algorithm  was  the  most  efficient  at  that  time.  Other 
shortest  path  routines  have  been  developed  since  1969,  e . g .  Yen  [160] 
and  O' Regan  [122],  but  while  efficient,  they  are  not  general.  Moreover 
these  algorithms  are  concerned  with  the  shortest  path  through  a  directed, 
cyclic  network. 

The  network  developed  for  this  research  is  directed  but 
contains  no  loops.  The  above-mentioned  algorithms  could  be  used  in 
this  application,  but  an  easier  approach  exists  since  the  network  can 
be  numbered  to  satisfy  the  I~J  rule.  Since  only  those  activities  whose 
initial  node  equals  or  exceeds  the  origin  node  of  the  path  need  to  be 
evaluated,  half  the.  number  of  additions  and  comparisons  is  required  in 
an  acyclic  network  compared  to  a  cyclic. 

By  identifying  the  optimal  predecessor  node,  the  initial  node 
producing  the  lowest  duration,  the  optimal  path  can  be  found  by  working 
backwards  through  the  network. 

Elmaghraby  [42,  p.  10]  proposed  a  dynamic  programming  algorithm 
for  the  shortest  path  through  a  directed,  acyclic  network  which  uses  the 
following  relationship: 

a.  =  min  .  .  (a.  +  D..)  (20) 

J  x  ->  j  l  xq 

where:  a.  =  the  duration  of  the  shortest  path  at  node  j 

D. .  =  the  duration  of  the  path  between  nodes  i  and  j. 
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The  following  are  the  main  steps  in  the  formulation  of 

Elmaghraby’s  algorithm  used  in  LOGNAP : 

STEP  (a)  -  identify  the  minimum  possible  terminal  number  node  for  an 
activity.  Note  that  this  is  the  initial  node  plus  one, 
since  the  terminal  node  of  an  activity  must  exceed  the 
initial  node  number  to  satisfy  the  I-J  rule . 

STEP  (b)  -  set  the  optimum  predecessor  to  zero  and  the  path  cost  to 
infinity. 

STEP  (c)  -  consider  each  activity  to  see  if  its  terminal  node  is  the 
minimum  possible,  and  if  its  initial  node  of  the  activity' 
is  not  less  than  the  initial  node  for  this  path.  If  not, 
continue  Step  (c)  for  the  remaining  activities. 

STEP  (d)  -  test  if  the  path  length  from  the  initial  node  to  this  node 
is  minimal.  If  so,  establish  the  new  minimal  path  length 
and  optimal  predecessor.  Continue  at  Step  (c)  with  the 
remaining  activities. 

STEP  (e)  -  when  all  activities  have  been  considered,  establish  a  new 
initial  node,  and  continue  at  Step  (a) .  When  all  nodes 
greater  than  the  initial  node  of  the  path  have  been 
examined,  the  length  of  the  shortest  path  is  known. 

STEP  (f)  -  the  nodes  on  the  shortest  path  are  found  by  working 
backwards  from  the  terminal  node  to  the  origin  using 


the  optimal  predecessors  of  each  node. 


. 
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4.  RESULTS  AND  DISCUSSION 

To  illustrate  the  application  of  shortest  path  network 
analysis  to  the  selection  of  least-cost  forest  harvesting  systems, 
several  LOGNAP  runs  were  made.  These  runs  investigated  the  optimal 
and  near-optimal  harvesting  systems  under  a  set  of  standard  environ¬ 
mental  conditions.  The  effects  that  changes  in  these  conditions  had 
on  the  systems  selected  were  also  studied. 

The  expanded  network  diagram  shown  in  Figure  7  was  used  for 
the  runs.  The  network  consisted  of  133  activities  and  77  nodes,  and 
had  848  alternative  paths.  The  basic  data  for  the  activities  comprising 
the  network  including  initial  and  final  node  numbers,  activity  descriptions 
and  various  costing  variables  outlined  in  Section  3.2  are  summarized  in 
Appendix  D. 

The  standard  environmental  conditions  common  to  all  activities 
were  separated  into  economic  and  operating  variables.  These  standard 
parameters  are  summarized  in  Table  7.  Except  where  noted  in  the 
following  sections  these  values  for  the  variables  were  used. 

The  sections  which  follow  describe  the  results  of  the  runs 
and  suggestions  for  future  work. 
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Table  7.  Values  of  Variables  Describing  Standard  Environmental 
Conditions  Used  in  the  Sample  Analysis. 


ECONOMIC  VARIABLES 

Interest  Rate  (%)  . . .  12 

Fuel  Price  ($/litre)  .  0.22 

Repair  Time  In-Shift  (%  of  total  repair  time)....  85 

Fringe  Benefit  Factor  (%  of  hourly  wage)  ........  35 

Insurance  Rate  Factor  (%  of  unrecovered  value)  ..  3 
Rate  of  Lost  Utility  (%/year)  .  30 

OPERATING  VARIABLES 

Tree  Species  . . . .  Black  Spruce 

Average  Tree  Size  (DBH  in  cm)  .  18 

Stand  Density  (m3/hectare)  . .  100 

Average  Primary  Transport  Distance  (m)  . 150 

Hauling  Distance  -  Road  Class  1  (km)  . 96 

Hauling  Distance  -  Road  Class  2  (km)  . 16 

Hauling  Distance  -  Road  Class  3  (km)  .  3 

Average  Slope  (%) 
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4.1  Optimal  and  Near-Optimal  Harvesting  Systems 

LOGNAP  was  used  to  select  the  50  least-cost  forest  harvesting 
system  from  the  848  possible  alternatives.  Appendix  E.l  contains  the 
complete  computer  output  for  this  run. 

Fifty  paths  were  selected  so  as  to  minimize  computing  time 
and  storage  requirements.  It  was  felt  that  identifying  additional 
systems  beyond  this  number  would  provide  little  benefit  other  than 
satisfying  curiosity.  If  needed,  the  costs  of  non-selected  systems 
can  be  calculated  using  the  activity  duration  data  from  LOGNAP’ s  output. 

On  the  Hewlett-Packard  9845B  mini-computer  used  for  this 
research,  a  run  time  of  approximately  two  hours  was  required  to  trace 
the  50  paths.  This  relatively  long  computer  run  time  is  attributable 
to  the  small  size  of  the  computer  rather  than  algorithmic  inef f iciences . 
A  test  run  of  LOGNAP  on  a  larger  computer,  an  H.-P.  1000,  suggests  that 
run  times  several  hundred  times  shorter  can  be  attained  with  the 
programs. 

Results  of  this  run,  specifically  the  range  in  woods  costs 
within  the  top  50  systems  and  the  actual  systems  selected,  are  discussed 
in  the  next  two  sections. 

4.1.1  Wood  Costs 

Wood  cost  data  for  the  first  50  least-cost  alternative 
harvesting  systems  are  plotted  in  Figure  9.  The  systems  are  ranked 
in  ascending  order  based  on  wood  costs  measured  in  dollars  per  cubic 
metre  of  tree-length  wood  unloaded  at  the  mill  site.  Note  that  these 
costs  are  based  solely  on  machine  rates  and  do  not  include  such  costs 


as  roads,  supervision  and  overheads. 
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The  cost  curve  in  Figure  9  indicates  that  the  two  systems 
which  produced  the  lowest  wood  costs  are  significantly  cheaper  than 
those  subsequently  selected  for  the  standard  conditions.  The  two 
best  systems  have  wood  costs  10  percent  lower,  about  $l/m3,  less 
than  the  third  best  system.  Perhaps  more  important  however  is  that 
the  range  in  wood  costs  between  System  10  and  System  50  is  also 
about  $l/m3.  A  difference  of  only  $0.40/m3  existed  in  the  wood 
costs  of  System  26  and  System  50. 

In  other  words,  under  the  conditions  modelled  in  this  run, 
wood  costs  are  relatively  insensitive  to  the  actual  system  used.  The 
majority  of  alternative  systems  are  not  significantly  cheaper  than  other 
systems.  Figure  10  shows  that  80%  of  the  alternative  systems  produce 
only  about  one  third  of  the  total  difference  in  woods  costs  between 
System  1  and  System  50. 

Under  these  circumstances,  the  ability  of  shortest  path 
network  analysis  to  identify  both  optimal  and  non-op timal  is  a  great 
advantage.  The  penalty  costs  of  operating  some  machines  outside  their 
optimum  range  appear  to  be  so  small  that  a  company  might  be  prepared  to 
absorb  them  rather  than  make  additional  capital  expenditures.  As  such 
identifying  the  ranking  and  costs  of  sub-optimal  systems  may  be  of  equal 
importance  to  finding  the  least-cost  system. 
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Figure  10,  Frequency  distribution  of  wood  costs  among  the  50  least-cost  alternative  systems. 
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4.1.2  Systems  Selected 

Table  8  summarizes  the  alternative  subsystems  among  the 
50  least-cost  systems  as  a  means  of  identifying  the  types  of  equipment 
which  offer  the  best  alternatives.  Total  frequency  of  selection  for 
each  subsystem  are  located  at  the  bottom  of  the  table.  Note  that  the 
subsystems  which  occur  between  landing  and  mill  have  not  been  included 
since  they  are  common  to  all  systems. 

The  majority  of  selected  systems  involve  the  subsystems  of 
felling  and  bunching,  primary  transport  of  bunched  full-trees,  and 
delimbing  at  roadside.  No  particular  piece  of  equipment  appears  to 
be  superior  for  fulfilling  the  fell  and  bunch  function.  Grapple  and 
clam-bunk  skidders  appeared  to  be  superior  for  the  primary  transport 
stage.  Roadside  delimbing  with  flails  was  clearly  superior  but  the 
model  does  not  include  any  penalty  costs  for  poor  delimbing  quality. 

The  first  five  least-cost  systems  utilize  a  "swathing” 
machine  presently  under  development  which  fells  and  primary  transports 
combined  with  roadside  flail  delimbing.  This  machine  offers  signifi¬ 
cantly  lower  wood  costs  than  presently  available  systems  under  the 
conditions  where  the  machine  can  operate.  Somewhat  surprising  was 
the  low  ranking  of  systems  using  the  similar  feller-forwarder  concept 
(34th)  although  the  machine  has  recently  enjoyed  considerable  commercial 


success . 
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Table  8.  Main  Subsystems  in  the  50  Least-Cost  Alternative  Harvesting  Systems 
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4.2  Sensitivity  Analysis 

The  operating  environments  of  forest  harvesting  systems  are 
extremely  variable.  Economic  conditions  are  continually  changing.  The 
physical  environment  is  seldom  constant  even  within  a  cut-block.  Other 
conditions  set  by  management  decisions  are  subject  to  change. 

Given  this  variable  operating  environment,  identifying  the 
optimal  and  near-optimal  systems  is  of  limited  value  since  the  model 
parameters  are  fixed  for  any  run.  These  parameters  may  be  the  best 
estimates  of  present,  and  future  operating  conditions,  but  they  remain 
estimates  with  uncertainty  associated  with  them.  Thus  the  solution  of 
a  practical  problem  is  incomplete  without  investigating  how  changes  in 
the  operating  environment  affect  the  optimal  solution. 

Sensitivity  analysis,  the  study  of  how  the  parameters  affect 
the  optimal  solution  can  be  undertaken  using  LOGNAP .  To  illustrate  how 
sensitivity  analysis  could  be  applied,  several  potentially  important 
parameters,  energy  prices  and  tree  size,  were  studied.  Sensitivity 
analysis  was  also  used  to  determine  the  optimal  operating  zones  for 
harvesting  systems. 

4.2.1  Energy  Prices 

The  effects  of  energy  prices  were  chosen  for  evaluation 
because  of  their  volatile  nature  and  potential  scarcity.  An  analysis 
of  the  historical  impact  of  energy  expenditures  on  wood  costs  suggested 
that  wood  costs  increase  at  about  half  the  rate  of  energy  price  hikes 
on  a  constant  dollar  basis.  This  figure  was  used  to  account  for  the 


impact  of  different  fuel  prices  on  all  components  of  machine  costs. 
Fuel  prices  were  then  varied  from  $0. 165/litre  and  $0. 44/litre,  the 
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upper  figure  being  approximately  double  present  costs  while  holding 
the  other  standard  variables  constant.  Results  of  these  LOGNAP  runs 
are  summarized  in  Appendix  E.2. 

The  selection  sequence  of  the  standard  systems  was  greatly 
affected  by  changes  in  fuel  prices.  Table  9  compares  the  selection 
sequence  of  alternative  harvesting  systems  at  various  levels  of  fuel 
costs.  The  systems  are  numbered  based  on  their  ranking  under  the 
standard  conditions. 

While  the  three  least-cost  systems  maintained  their  rankings 
over  the  range  of  prices  studied,  some  systems  proved  to  be  very 
sensitive  to  energy  price  changes.  For  example,  a  manual,  "cut  and 
skid"  system  which  ranked  19th  under  standard  conditions  was  the  fourth 
cheapest  system  when  fuel  prices  were  doubled.  A  system  that  was 
unranked  when  fuel  cost  $0. 22/litre  was  selected  eighth  when  fuels 
cost  $0. 44/litre.  Other  so-called  manual  systems  showed  similar 
dramatic  improvements  in  their  rankings. 

The  overall  effect  of  a  doubling  in  energy  prices  was  a  67% 
increase  in  wood  costs  for  the  least-cost  system.  The  overall  shape 
of  the  system  cost  curve  (Figure  10)  remained  similar  but  was  shifted 
upwards  with  higher  energy  prices. 

Undertaking  sensitivity  analysis  with  LOGNAP  was  easy. 

Only  one  equation  and  one  standard  variable  required  changes  to 
determine  the  effects  of  energy  prices  on  the  selection  of  the 


least-cost  systems. 
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Table  9.  Changes  in  the  Selection  Sequence  of  the  Standard  Systems  under  Different 
Fuel  Prices 
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4.2.2  Average  Tree  Size 

One  of  the  most  important  operating  variables  affecting 
system  performance  is  average  tree  size.  It  enters  activity  duration 
calculations  in  the  cycle  time  equations  for  many  machines  and  in  the 
hourly  productivity  equation  (#18) .  To  study  the  sensitivity  of 
systems  selection  to  this  variable,  a  series  of  LOGNAP  runs  were  made 
varying  tree  diameters  at  breast  height  between  12  and  30  centimeters. 
Results  of  these  runs  are  summarized  in  Appendix  E.3. 

Figure  11  illustrates  the  observed  effects  of  average  tree 
size  on  the  least-cost  system  selection.  While  the  swather  with  roadside 
flail  delimbing  remained  optimal  over  most  of  the  range  in  tree  size, 
a  tractor-mounted  shear  system  was  the  least-cost  s3/stem  with  trees 
greater  than  26  cm  DBH.  This  system  however,  exhibited  great  sensitivity 
to  average  tree  size. 

A  swather-based  system  with  a  roll  delimber  was  optimal 
with  very  small  trees  though  very  marginally  ($0.01/m3).  In  actual 
practice,  it  is  unlikely  that  this  system  would  be  utilized  since  a 
stand  of  trees  with  a  12  cm  diameter  is  on  the  borderline  of 
merchantability. 

An  advantage  of  the  shortest  path  network  analysis  approach 
to  system  selection  is  apparant  from  this  example  of  sensitivity 
analysis.  Prediction  equations  of  non-linear  form  are  easily  incor¬ 
porated  into  the  model  without  the  need  for  approximation  using  line 
segments  as  mathematical  programming  approaches  would  require.  Three 
or  four  segments  would  be  necessary  to  accurately  describe  the  cost 
curve  of  the  tractor-mounted  shear  system,  but  only  one  equation  was 


needed  using  LOGNAP. 
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WOOD  COST  ( $/m3 ) 


Figure  11. 


AVERAGE  TREE  SIZE 
(DBH  in  cm) 


Effects  of  average  tree  size  on  harvesting  costs  of 
the  least-cost  harvesting  system. 
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4*2.3  Optimal  Operating  Zones 


One  potential  application  of  sensitivity  analysis  with  LOGNAP 
is  to  identify  the  optimal  operating  zones  for  harvesting  systems.  To 
illustrate  such  a  use,  sensitivity  analyses  were  conducted  on  a  range 
O"  average  tree  sizes  and  stand  densities.  These  parameters  were  varied 
across  the  range  of  conditions  typically  found  in  Eastern  Canada:  tree 
diameters  at  breast  height  from  12  to  30  centimeters  and  stand  volumes 
from  60  to  180  cubic  metres/hectare.  The  results  of  these  LOGNAP  runs 
are  summarized  in  Appendix  E . 3 . 

Both  stand  parameters  had  a  significant  impact  on  the  sequence 
of  system  selection.  Many  systems  which  had  not  been  selected  in  the 
previous  analyses  became  optimal  or  near-optimal.  Over  the  range  of 
parameters  examined,  wood  costs  more  than  doubled.  Table  10  summarizes 
the  wood  costs  for  the  three  systems  which  became  optimal  over  conditions 
studied. 

Three  systems  were  optimal  for  practical  purposes  under  some 
combination  of  stand  density  and  tree  size.  A  fourth  system  was  optimal 
with  very  small  trees  as  discussed  in  the  previous  section  but  will  not 
be  considered.  The  swather-roadside  flail  system  was  the  least-cost 
system  over  much  of  the  range.  A  tractor-mounted  shear  with  tracked 
cable  skidding  and  roadside  flail  delinking  was  optimal  with  large  trees 
and  low  densities.  A  manual  fell,  tracked  cable  skidder  and  roadside  flai 
system  was  the  least-cost  system  with  an  average  tree  dicimeter  of  18  cm 
and  a  stand  density  of  60  m  Vha.  Figure  12  combines  the  response 
surfaces  of  the  three  systems.  The  shaded  areas  outline  the  conditions 
under  which  each  system  is  optimal. 
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Table  10.  Wood  Costs  ($/m3)  for  the  Least-Cost  Harvesting  Systems 
Under  Varying  Stand  Volumes  and  Average  Tree  Sizes. 
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_j  Swather  System 
□  Tractor-Mounted  Shear  System 


Figure  12.  Least-cost  harvesting  systems  under  varying  stand  volumes 
and  average  tree  sizes. 
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4,3  Program  Modifications  and  Extensions 


Several  modifications  could  be  made  to  LOGNAP  to  improve  the 
accuracy  and  usefulness  of  the  program.  The  following  are  some  areas 
for  further  work. 

4.3.1  Costing  Models 

The  costing  model  used  in  this  research  while  more  realistic 
than  many  used  in  the  past  e . g .  [15]  could  still  be  improved.  The 
treatment  of  equipment  salvage  values  and  repair  costs  in  particular 
need  refinement. 

The  assumption  that  a  machine's  salvage  value  equals  what  is 
in  effect  its  undepreciated  worth  is  not  realistic  in  the  present  used- 
machinery  market.  Price  escalations,  largely  due  to  a  high  inflation 
rate,  have  recently  been  so  rapid  that  salvage  values  of  equipment 
trade-ins  may  equal  or  exceed  the  original  purchase  price .  Studies  of 
logging  equipment  salvage  values  over  time  are  not  available . 

A  similar  situation  exists  with  repair  costs.  Though  some 
companies  can  establish  long-term  repair  costs,  these  are  not  generally 
available.  This  research  assumed  a  long-term  average  cost  for  repairs. 
While  this  approach  gives  reasonable  estimates  of  total  expenditures , 
it  may  not  provide  a  good  picture,  of  the  actual  timing  of  these  costs . 

With  a  cash  flow- type  model,  timing  is  a  critical  concern. 

The  development  of  an  after-tax  cash  model  should  be  considered, 
although  machine  costing  is  presently  done  on  a  before-t:ax  basis  in  the 
Canadian  logging  industry.  This  approach  would  provide  a  better  approxi¬ 
mation  of  the  company's  actual  cash  flow. 
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The  programming  of  several  alternative  costing  models  so  that 
users  could  select  the  approach  they  wish  to  use  may  be  desireable. 

Cash  flow  models  might  also  incorporate  estimates  of  future  costs  for 
things  such  as  energy  and  inflation. 

4.3.2  Productivity  Equations 

While  the  equations  used  to  predict  productivity  represent  the 
best  available  information,  they  were  frequently  found  to  be  deficient. 
In  most  cases  for  example,  no  relationship  had  been  established  for  slop 
or  soil  conditions  though  such  must  exist  and  exercise  a  major  influence 
on  the  system  selected.  Most  machine  studies  do  not  examine  sufficient 
stand  and  terrain  conditions  to  quantify  the  impact  of  these  variables 
on  cycle  times . 

In  other  cases  the  equations  used  were  taken  from  studies  of 
prototype  machines  and  old  models.  Whether  these  relationships  apply  to 
present  equipment  is  debatable.  Certainly  the  establishment  of  accurate 
relationships  between  cycle  times  and  various  operating  conditions  is 
necessary  to  minimize  the  selection  of  ncn-optimal  systems. 

An  alternative  would  be  to  model  the  physical  relationships 
involved  with  a  working  machine  and  derive  cycle  time  estimates  by 
calculation  rather  than  from  observation. 

4.3.3  Data  Set  Edit  Routine 

No  changes  were  made  in  the  basic  data  set  during  the  sample 
application.  As  presently  programmed,  such  changes  could  only  be  made 
by  complete  re-entry  of  the  data  using  the  INPUT  program.  A  routine  for 
editing  the  data  set  would  be  a  useful  addition  to  LOGNAP .  This  would 
permit  the  updating  of  purchase  prices  and  labour  rates ,  and  evaluation 
of  the  impact  of  machine  specific  variables  on  systems  selection. 
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4.3.4  Uncertain  Activity  Durations 

To  compensate  for  imprecise  costing  and  productivity  estimates, 
uncertainty  could  be  incorporated  into  the  estimates  of  activity  duration. 
A  probability  distribution  for  the  estimated  cost  of  each  activity  using 
optimistic,  likely  and  pessimistic  estimates  of  duration  as  is  done 
with  PERT  could  be  established.  The  expected  cost  and  a  variance  measure 
for  any  system  could  then  be  obtained.  If  a  large  variance  is  associated 
with  the  shortest  path,  or  only  small  differences  exist  in  the  length  of 
the  shortest  and  sub-shortest  paths,  improved  system  selection  might  be 
achieved  by  considering  risk  and  minimizing  the  uncertainty. 

4.3.5  Negative  Activity  Durations 

One  advantage  of  Yen's  algorithm  is  that  the  concept  of  a 
negative  activity  duration  can  be  handled.  This  would  permit  compensa¬ 
tion  for  non-equivalent  activities  to  be  included  within  the  network. 

For  example  shear  damage  to  tree  butts  may  produce  a 
significant  loss  of  value  (10-15%)  for  fibre  destined  for  lumber  [50], 
but  little  loss  in  value  if  intended  for  pulping.  Penalizing  shearing 
equipment  using  a  negative  duration  would  place  both  the  equipment  and 
product  on  an  equivalent  base.  Similarly  an  expanded  network  of  all 
types  of  harvesting  systems  could  be  prepared  with  a  compensating 
activities  to  place  all  systems  to  a  common  basis.  Limbing  quality 
is  another  area  where  the  concept  of  penalty  costs  could  be  incorporated. 
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4.3.6  Optimal  Operating  Layouts 

Comparisons  made  with  LOGNAP  at  present  are  based  on  all 
systems  operating  under  the  same  logging  layouts.  While  this  may  be 
necessary  for  selecting  systems  for  particular  situations,  costs 
when  all  equipment  is  operating  under  optimal  layouts  can  be  useful 
particularly  for  planning  purposes. 

The  treatment  of  average  primary  transport  distance  in  the 
present  model  illustrates  the  problem.  A  constant  distance  was  assumed 
for  all  primary  transport  equipment.  However  if  this  distance  is  not 
optimal  for  a  particular  machine,  any  system  using  that  machine  could 
produce  wood  at  a  lower  cost  with  appropriate  layout  changes.  At  the 
same  time,  road  costs  could  be  incorporated  in  the  model. 

Modification  of  the  program  to  permit  selection  of  systems 
operating  under  optimal  layouts  is  desirable. 

4.3.7  System  Balancing 

LOGNAP  identifies  the  equipment  which  would  produce  the  lowest 
cost  wood.  It  does  not  consider  the  problem  of  how  many  of  each  machine 
type  is  necessary  to  meet  organizational  goals.  Further  work  could 
incorporate  the  optimal  system  in  a  mixed-integer  mathematical  program¬ 
ming  model  which  balances  the  system  within  overall  company  objectives. 
Several  forest  harvesting  simulators  which  are  capable  of  modeling  most 
systems  e. g.  [69;  152]  could  also  be  utilized  to  balance  the  system. 
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5.  SUMMARY  AND  CONCLUSIONS 


Alternative  methods  of  delivering  tree-length  wood  to  a 
processing  plant  were  outlined  as  a  directed,  acrylic  network  diagram. 
Completed  operations  were  represented  by  nodes  in  the  network.  Arcs 
joining  the  nodes  described  possible  sequences  of  activities  for 
fulfilling  the  system’s  function. 

The  duration  of  each  activity,  measured  in  dollars  per  cubic 
metre  of  wood  produced,  was  estimated  using  the  hourly  costs  and 
productivity  of  the  equipment  and  labour  involved.  Hourly  costs  were 
calculated  on  an  annual  equivalent  basis  using  a  before-tax  cash  flow 
model  of  capital  and  operating  costs.  Productivity  estimates  used 
cycle-time  prediction  equations  derived  from  time  studies,  and 
parameters  such  as  tree  and  load  sizes  to  characterize  the  operating 
environment  and  the  equipment. 

A  package  of  computer  programs  was  developed  to  select  the 
optimal  (least-cost)  and  near-optimal  systems  using  shortest  path 
network  analysis.  The  systems  were  identified  using  Yen’s  algorithm 
[159 J  to  find  a  specified  number  of  shortest  paths.  The  shortest 
distance  and  nodes  on  a  path  between  some  node  and  the  network’s 
terminal  node  were  found  using  a  dynamic  programming  algorithm 
proposed  by  Elmaghraby  [42]. 

Of  the  alternatives  modelled,  those  systems  using  a  feller- 
forwarder  based  on  the  swathing  principle  produced  the  lowest  wood 
cost  under  most  conditions.  A  system  consisting  of  a  swather  and 
flail  delimbing  at  roadside  was  the  least-cost  system  over  a  wide 
range  of  energy  prices,  tree  sizes  and  stand  volumes.  A  tractor- 
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mounted  shear  system  produced  the  lowest  cost  wood  under  conditions  of 
low  stand  volume  and  large  tree  size.  The  widely  used  manual  felling 
system  was  the  least-cost  system  with  trees  diameters  of  about  24  cm 
and  low  stand  volumes. 

Shortest  path  network  analysis  proved  to  be  an  excellent  means 
of  identifying  least-cost  forest  harvesting  systems.  A  network  diagram 
in  itself  is  a  useful  aid  to  visualizing  the  sequence  of  activities  and 
overlooked  alternatives  for  meeting  a  system’s  function.  Given  the  small 
differences  observed  between  the  wood  costs  of  most  systems,  the  identifi¬ 
cation  of  more  than  the  least-cost  system  is  highly  desirable.  Sensitivity 
analysis  is  easily  and  quickly  completed. 

Perhaps  the  greatest  advantage  of  shortest  path  network  analysis 
is  that  a  great  number  of  alternatives  can  be  defined  in  terms  of  standard 
estimating  equations,  minimizing  the  time  spent  defining  the  problem  and 
maximizing  the  understanding  of  non- technical  personnel. 
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APPENDIX  A 

PRICE  INDICES  USED  TO  DEFLATE  COST  STATISTICS 

A. 1  Consumer  Price  Index 

The  Consumer  Price  Index  (CPI)  is  applied  to  total  annual 
wages,  so  that  employee  earnings  can  be  compared  on  an  equal  buying- 
power  basis  over  the  time  frame.  CPI  values  for  1960  to  1978  are 
tabulated  in  Table  A.l.  These  values  were  applied  to  calculate  the 
"Average  Wage"  column  in  Table  1. 

A. 2  Forestry  Machinery  and  Equipment  Index 

The  Forestry  Machinery  and  Equipment  Index  (FMEI )  is  applied 
to  capital  expenditures  on  logging  activities  to  permit  the  comparison 
ot  capital  investments  made  at  different  times  on  an  equal  basis.  FMEI 
values  are  tabulated  in  Table  A.l.  These  values  were  used  to  equate  the 
capital  expenditures  in  Figure  4. 

A. 3  Implicit  Price  Indices  for  Gross  National  Product 

The  Implicit  Price  Indices  (IPI)  permit  the  comparison  of 
costs  composed  of  many  indirectly  related  items.  The  effects  of  all 
variables  included  in  the  Gross  National  Product  are  averaged,  permit¬ 
ting  unrelated  costs  to  be  compared  on  a  common  basis.  IPI  values  are 
tabulated  in  Table  A.l.  These  values  were  applied  to  standardize  wood 
costs  in  Figure  3,  and  the  productivity  indices  in  Table  2. 


. 


' 


Table  A.l.  Price  Indices  for  the  Years  1960  to  1978  [146;  148;  151]. 


YEAR 

CONSUMER 

PRICE 

INDEX 

FORESTRY 
MACHINERY  AND 
EQUIPMENT 
INDEX 

IMPLICIT 

PRICE 

INDEX 

1960 

74.3 

83.6 

71.0 

1961 

75.0 

84.8 

72.4 

1962 

75.9 

87.3 

73.4 

1963 

77.2 

90.4 

74.8 

1964 

78.6 

96.4 

76.6 

1965 

80.5 

100.7 

79.1 

1966 

83.5 

101 . 0 

82.6 

1967 

86.5 

94.6 

85.9 

1968 

90.0 

91.5 

88 . 7 

1969 

94.1 

94.2 

92.6 

1970 

97.2 

97.4 

96.9 

1971 

100.0 

100.0 

100.0 

1972 

104..  8 

102.8 

104.8 

1973 

112.7 

107.1 

113.4 

1974 

125.0 

122.2 

129.4 

1975 

138.5 

142.2 

148.8 

1976 

148.9 

149.0 

163.3 

1977 

160.8 

161.7 

171.4 

1978 

175.2 

180.7 

182.3 
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APPENDIX  B 

EQUATIONS  USED  TO  CALCULATE  CYCLE  TIMES 


B.l  Formulae  for  Common  Variables 

Percentage  of  grapple  area  occupied  by  full-tree  stems  [82]. 

PACKINGFT  (decimal  %)  =  0.188  +  0.0213  x  DSH  (cm) 

Percentage  of  grapple  area  occupied  by  tree-length  stems  [82 J. 

PACKINGTL  (decimal  %)  =  0.440  +  0.0213  x  DSH  (cm) 

Average  travel  distance  (DIST) 

DIST  (m)  -  Wander  Factor  x  Average  Primary  Transport  Distance  (m) 
=  1.4  x  APTD 


Number  of  trees  per  hectare  (NTPH) 


NTPH 


Merchantable  Volume  per  Hectare  (m3/ha) 
Merchantable  Volume  per  Tree  (m^ ) 


Machine  payload  [154;  pp.  60-64] 


Payload 


T  (TE  -  R  x  TW) 

FS  -  T  x  R  x  P 


where:  T  =  traction  coefficient 

TE  =  tractive  effort  (kg) 

FS  =  skidding  coefficient 

R  =  rolling  resistance  coefficient 
TW  =  tare  weight  (kg) 

P  =  proportion  of  payload  supported  by  skidder 

Number  of  trees  per  cycle 

NT  _  Payload  (kg) _ 

Weight  per  Tree  (kg) 


or 


, 

{  .  iiH  -I J  |  '  _ 


■ 


Integer  (Rated  Boom  Capacity  at  Maximum  Reach  (kg). 

3  x  Weight  per  Tree  (kg) 

or 


Integer  f4  x  Grapple  Area  (cm2)  x  PACKING 
6  tt  x  DBHZ  (cm) 


B.2  Equations  for  Various  Tree  Parameters 

Regression  equations  for  the  tree  parameters  needed  in  the 
cycle  time  prediction  equations  are  contained  in  Tables  B.l  through 
3.5  for  the  four  main  softwood  species  of  Eastern  Canada, 


Regression  Equations  for  Calculating  Stump  Diameter  (15  cm  height)  [2,  p.  8-11]. 
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B.3  Cycle  Time  Equations 

1.  Feller-buncher ,  wheeled,  shear  -  Forano  BJ-20  [131] 

CT (1)  =  50  +  212  x  MV 

2.  Feller-buncher,  tracked,  shear  -  Drott  40  [127] 


3, 


4, 


CT(2)  =  40  +  21.2  x  MV 


Feller-buncher,  multiple-tree  shear  -  Koehring  KFB4  [85] 

TIBS 

NT  =  INT  (  — )  +  1 

v3  x  WTFT; 

CT(3)  =  49.4  -  7.2  x  NT  +  42.4  x  MV 
Feller-buncher  wheeled,  shear  -  Melroe  Bobcat  1075 


CT (4)  —  43  + 


3460 

NTPH 


5.  Feller-buncher,  wheeled,  saw  -  Kockums  880  [10;  46] 

CT (5)  =  40  x  e  (DBH/1065) 

6.  Feller-buncher,  tracked,  saw  -  Drott  40  [50;  127] 

CT (6)  =  40  +  113  x  MV 

7.  Feller-buncher,  tracked,  auger  -  Drott  40  [51;  127] 

CT(7)  =  40  +  113  x  MV 

8.  Feller-buncher,  wheeled,  multiple-tree  saw  -  Hydro-Ax  Swather 

1600 


NT  =  INT  (- 


+  1 


'3  x  WTFT' 

CT (8)  =  70  +  0.70  x  SLOPE  -  7  x  NT 

9.-16.  Dummies 

17.  Fell  -  man  with  chain  saw  [154] 

CT (17 )  =  0.650  x  DBH  x*81 

18.  Fell,  tracked  shear  -  Caterpillar  D6  [29] 

pths'i  -  _ 600u _ 

V  ;  95.46  +  0.085  x  NTPH 

19.  Fell,  tracked  saw  -  FMC  200  with  feller-director  [84] 


CT (19)  =  60  +  116.5  x  MV 


■ 

. 


. 


. 


Ill 


20. 


Fell,  tracked  saw  -  Caterpillar  D6  [29] 


CT(20) 


6000 _ 

90  +  0.085  x  MV 


21.-24.  Dummies 

25.  Delimb  stumpside  (loose)  -  Logma  T-310  [129] 

CT(25)  =  27  +  ||||3  +  24.7  x  MV 

26.  Delimb  stumpside  (loose)  -  Drott  40  with  roll  delimber 

CT (26)  =  53  +  14  x  MV 


[45] 


27.-29.  Dummies 

30.  Fell  and  delimb  -  man  with  chain  saw  [154] 

CT (30)  =  0.55  x  DBH  1,81  +  0.17  x  DBH  2,42 

31.  Fell  and  delimb  -  John  Deere  743  [49] 

CT(31)  =  58  +  0.86  x  ML 


32.-33.  Dummies 

34.  Fell,  delimb  and  bunch  -  Timber jack  TJ-30  [130] 

CT (34)  =  66  -  0.60  x  SLOPE  +  102.4  x  MV 

35.  Fell,  delimb  and  bunch  -  Koehring  KFB4  with  Timmins  head  [85;  45] 

CT (35)  =  66.4  +  14  x  MV 


36.  Fell,  delimb  and  bunch  -  Drott  40  with  Timmins  head  [45;  47] 

CT(36)  =  62  +  98  x  MV 

37. -39.  Dummies 

40.  Delimb  stumpside  (bunched)  -  Logma  T-310  [129] 

CT (40)  =  35  +  102.4  x  MV 

41.  Delimb  stumpside  (bunched)  -  GLFP/NESCO  flail  [66] 

CT(4i)  =  40 


42.  Delimb  stumpside  (bunched)  -  Drott  40  with  Timmins  head  [45] 
CT (42)  =  30  +  38.5  x  MV 


43.-45.  Dummies 


. 


V 

■ 


■ 


112 


46.  Primary  transport  (bunched  full  trees),  wheeled,  cable  skidder 
John  Deere  540  [56] 


47. 


48. 


0.55  (6375  -  0.50  x  7580)  _ 

0.65  -  0.55  x  0.5  x  0.5  "  l'11^  kg 


Payload 
NT 

CT(46)  =  32  +  229  x  MV  + 


2775 
INT  (——--) 
VWTFT; 


200  +  2.5  x  DIST 
’  NT 


D L -ty  transport,  (ounched  full  trees)  ,  wheeled,  cable  skidder  — 
John  Deere  740  [56] 

Payload  =  J113g5  ~  =  5685  kg 

0.65-0.55x0.5x0.5  b 

NT  =  INT  (|£§§) 

CT(47)  =  32  +  229  x  MV  + 

NT 

Primary  transport  (bunched  full-trees),  wheeled,  grapple  skidder 
John  Deere  540  [19] 


49. 


NT 


_  TNT  (-4  x  8082  x  PACKINGFT, 


CT(48)  =  22.6  x  MV  + 


it  x  DBlI7 

322  +  2.5  x  DIST 
NT 


Primary  transport  (bunched  full-trees),  wheeled,  grapple  skidder 
John  Deere  740  [19] 


INT  (4  x  13935  x  PACKINGFT^ 


NT 

CT(49)  =  22.6  x  MV  + 


7T  x  DBH^ 

400  +  2.5  x  DIST 
NT 


50. 


51. 


Primary  Transport  (bunched  full-trees),  wheeled,  clam  skidder 
Timber jack  520  [128] 


TNT  (4  x  32515  x  PACKING?!, 


NT 

CT (50)  =  162  x  MV  + 


v  x  DBH2  ' 

328  +  2.85  _x_ DIST 

NT 


Primary  Transport  (bunched  full-trees),  wheeled,  forwarder 
Koehring  KF2  [85;  48;  65] 


NT  =  INT  (“§§--) 


CT (51)  =  35  + 


WTFT 

52  +  5,84  x  DIST 
NT 


. 

r 


< 


. 


■ 


■ 


113 


52.  Primary  Transport  (bunched  full-trees),  tracked,  cable  skidder 
FMC  200  [84;  132] 


NT  =  INT 

wift' 


CT (52)  = 


1200  +  2.36  x  DIST 
NT 


53. 


Primary  Transport  (bunched  full- trees) ,  tracked,  clam  skidder 
FMC  200  [84:  132] 


NT 

CT(53) 


_  A  x  13285  x  PACKINGFTs 
1  x  DBH2  } 

_  1050  +  2.36  x  DIST 

NT 


54.-61.  Dummies 

62.  Primary  Transport  (loose  full-trees),  wheeled,  cable  skidder 
John  Deere  540  [56] 


63. 


27/S 

NT  =  INT  (f~) 

WTFT 

CT(62)  =  40  +  229  x  MV  + 


257  +  2.5  x  DIST 
NT 


Primary  Transport  (loose  full- trees) ,  wheeled,  cable  skidder 
John  Deere  740  [56] 


NT  -  INT 

V'WTFT' 

CT (63)  =  40  +  229  x  MV  +  -x  DI— - 

NT 

64.  Primary  Transport  (loose  full-trees),  wheeled,  clam  skidder 
Timber jack  520  [128] 

NT  =  int  (4  x  32515  x  PACKINGFT^ 


tt  x  DBH 

CT(64)  =  200  x  MV  +  600  +  2~85 

NT 


65.  Primary  Transport  (loose  full-trees),  wheeled  forwarder 
Koehring  KF2  [85] 


NT 


INT  (——---) 
VWTFT  ; 


CT(65)  -  55  +  -2L  -DI.ST 

v  ;  NT 


, 


- 


66. 


114 


Primary  Transport  (loose  full-trees),  tracked,  cable  skidder  ~ 
FMC  200  [84;  132] 


NT 


=  INT  (■ 


6200 

WTFT 


) 


CT (66) 


1350  +  2.36  x  DIST 
NT 


Primary  Transport  (loose  full- trees) ,  tracked,  clam  skidder  - 
FMC  200  [84;  132] 


NT 

CT (67) 


™  ,4  x  13285  x  PACKINGFT 
IN1  (  tt  x  DBH2 - 

1300  t  2.36  x  DIST 

NT 


68.-73.  Dummies 

74.  Fell  and  Primary  Transport,  shear,  skid  -  Timber jack  520  with 
head  L128;  45] 


NT  =  INT  (4  x  3251.5  x  PACKINGFT ^ 


CT(74)  =  50  + 


7T  x  DBHZ 

1000  +  2.85  x  DIST 
NT" 


75, 


Fell  and  Primary  Transport,  multiple-tree  shear,  forward 
Koehring  KFF  [85;  48;  65] 


NT  =  INT 

WTFT  7 

rrx'(nr\  /o  /  52  +  5.84  x  DIST  , 

Cl  (75)  =  49.4  +  - — -  +  42.4  x  MV 

1135 

-  7.2  x  (INT  (  --■ 1+1) 

3  x  WTFT7  7 

76.  Fell  and  Primary  Transport,  multiple- tree  saw,  forward 
Koehring  KFF  [85] 

NT  ~  TNT  ( 

NT  -  INT  ) 

n<vr-7£\  _  /n  /  .  52  +  5.84  X  DIST 

CT (76)  =  49.4  + - — - +  55  x  MV 

NT 


- 7-2  *  CINI  (3lr§f?T)+1) 


77. 


Fell  and  Primary  Transport,  saw,  forward  -  PAPCO  Swather 
CT (77)  -  10000/(1.02  x  NTPH) 


-81. 


Dummies 


Primary  iransport  (loose  tree— lengths) ,  x^heeled,  cable  slcidder 
John  Deere  540  [56] 


Payload 


0.55  (6375  -  0. 50  x  7580)  . 

0.55  -  0.55  x  0.50  x  0.~60_  ”  Jby°  kg 


NT 


INT  ( 


3690 

WTTL 


) 


CT(82) 


=  30  +  229  x  MV  + 


257  +  2.5  x  DIST 
NT . 


Primary  Transport  (loose  tree— lengths) ,  wheeled  cable  skidder  — 
John  Deere  740  [56] 


Payload  = 


0.55  (11365  -  0.50  x  12135) 

0.55  -  0.55  x  0.50  x  0.60 


7570  kg 


NT 


INT  (75Z0) 
wTTL ' 


CT(83)  =  30  +  229  x  MV  +  ‘A ■ 

NT 


Primary  Transport  (loose  tree— lengths) ,  wheeled  grapple  skidder 
John  Deere  540  [19] 


NT 


INT  /4  x  8082  x  PACKINGTL^ 


CT(84)  =  22.6  x  MV  + 


it  x  DB  H‘;- 

600  +  2.5  x  DIST 
NT 


Primary  Transport  (loose  tree-lengths),  wheeled,  clam  skidder  - 
Timber jack  520  [128] 

NT  =  tnt  f4  X  32515  x  PACKINGTL> 

’  ^  "VinDBH2’  '  "  ' 

CT (85)  =  175  x  MV  +  An85  x  DIST 


Primary  Transport  (loose  tree-lengths) ,  wheeled  forwarder  - 
Koehring  KF2  [85] 


=  INT  ( 


29000 

WTTL 


) 


CT(86) 


55 


52  +  5.84  x  DIST 

NT 


116 


Primary  Transport  (loose  tree-lengths),  tracked,  cable  skidder  - 
PMC  200  [84;  132] 


NT 


=  INT 


r6200 

'WTTL 


CT(87) 


1350  +  2,  36  x  DIST 
NT 


Primary  Transport  (loose  tree-lengths),  tracked,  clam  skidder  - 
PMC  200  [84;  132] 


NT 


=  INT  (- 


x  1.3285  x  PACKING!! 

IT  X  DBH2 


CT (88) 


1300  +  2.36  x  DIST 
NT 


89.-95.  Dummies 

96.  Primary  Transport  (bunched  tree-lengths),  wheeled,  cable  skidder 
John  Deere  540  [56] 


NT 

Cl (96) 


INT 

''WTTL 


30  +  229  x  MV  + 


200  a  2.5  x  DIST 
NT 


97.  Primary  Transport  (bunched  tree-lengths),  wheeled,  cable  skidder 
John  Deere  740  [56] 


NT 


CT(97)  =  30  a  229  x  MV  a 


200  a  2,5  x  DIST 
NT 


98.  Primary  Transport  (bunched  tree-lengths),  wheeled,  grapple  skidder 
John  Deere  540  [56] 


NT  =  int  (4  x  8082  x  PACKINGTL^ 


CT(98)  =  22.6  x  MV  + 


tt  x  DBHZ 

322  a  2.5  x  DIST 
NT 


99.  Primary  Transport  (bunched  tree-lengths),  wheeled,  grapple  skidder 
John  Deere  740  [56] 


NT  =  int  x  13935  x  PACKINGTL^ 


CT (99)  =  22.6  x  MV  + 


tt  x  DBHZ 

400  +  2.5  x  DIST 
NT 


■ 


■ 


' 


100.  Primary  Transport  (bunched  tree-lengths),  wheeled,  clam  skidder 
Timber jack  520  [128] 


NT 


,4  x  32515  x  PACKINGTL, 

1JN.L  { - — - - .) 


7T  x  DBH 


CT (100)  =  162  x  MV  + 


328  +  2.85  jx_  PI  ST 

NT 


101.  Primary  Transport  (bunched  tree-lengths),  wheeled  forwarder 
Ko eh ring  KF2  [85] 


INT  (— ~~) 


NT 

CT (101)  =  35  + 


WTTL 

52  +  5.84  x  DIST 
NT 


102.  Primary  Transport  (bunched  tree-lengths),  tracked,  cable  skidder 
FMC  200  [84;  132] 


NT 


-  TNT  f620CN 

-  INT 


CT(102)  = 


1200  +  2.36  x  DIST 
NT 


103.  Primary  Transport  (bunched  tree-lengths),  tracked,  clam  skidder 


FMC  200  [84;  132 


NT 

CT(103)  = 


,4  x  13285  x  PACKINGTL. 
=  INT  ( - — — ^552 - ) 


7T  x  DBH' 

1050  +  2.36  x  DIST 
NT 


104.-111.  Dummies 

112.  Delimb  roadside  -  Logma  T-310  [129] 

CT (112)  =  38  +  10.6  x  MV 

113.  Delimb  roadside  -  Harricana  [52] 

CT(113)  =  45 

114.  Delimb  roadside  -  Koehring  KBL 

CT(114)  =  50 

115.  Delimb  roadside  -  Drott  40  with  Timmins  roll  delimber  [47] 


CT (115)  -  20  +  49.4  x  MV 


. 


116.  Delimb  roadside  -  Hydro  Ax  500  flail  [53] 


CT(116)  =  18 

117.  Delimb  roadside  -  Hydro  Ax  flail  system  [53] 

CT (117 )  =  20 

118. -121.  Dummies 

122.  Loader  spreading  bunch  for  delimbing  -  Tanguay  14030  [53] 

CT (122)  =  20 

123.  Delimbing  touch-up  -  man  with  chain  saw 


CT (123) 

=  15 

124.  Harvest  and 

primary  transport  --  Koehring  RTL  [85] 

NT 

=  INI  ( 20000 ) 

VWTTL  ' 

CT (124) 

52  +  5.84  x  DIST 
=  /5  +  NT 

125.  Harvest  and  primary  transport  -  Koehring  KFF  with  Timmins  head 
[47;  85] 


NT 

= IKT 

CT(125) 

300  +  5.84  x  DIST 
-  54  +  49  x  MV  +  -  ■ 

INI 

126.-127.  Dummies 

128.  Load,  tracked  -  Caterpillar  235  [81] 


NT 

,4  x  1000  x  PACKINGTL, 

-  IN  I  (  •  "VP-TV72  ) 

ir  x  DBHZ 

CT (128) 

47.7 

NT 

129.  Dummy 

130.  Secondary  transport,  double  load  on  semi-trailer  [81] 


NT 

=  INI  (4100°) 

VWTTL 

CT (130) 

190  +  86  x  DISTRC1  +  144  x  DISTRC2  +  355  x  DISTRC3 

NT 

- 

V 


119 


131. 

132. 


Dummy 


Unload  -  Rago 
NT 

CT (132)  = 


Wagner  [81] 


INT 


41000 

WTTL 


300 

NT 


133.  Dummy 


APPENDIX  C 


COMPUTER  PROGRAMS 


C.l  Program  Listings 

The  programs  contained  in  this  appendix  make  up  LOG'NAP  . 
Internal  storage  limitations  of  the  Hewlett-Packard  9845.B  computer 
used  in  this  work  required  that  LOGNAP  be  broken  into  four  segments. 

The  first  program  called  "MAIN"  is  the  master  control 
program  controlling  the  branching  to  the  other  modules,  renumbering 
the  nodes  to  satisfy  the  I-J  rule,  and  outputting  the  final  shortest 
path  information. 

The  second  program  entitled  "INPUT”  is  used  to  input  and 
store  the  basic  network  description  and  costing  information  in  a 
file.  On  the  particular  installation  used  in  this  work,  data  was 
stored  on  a  floppy  disk. 

The  third  program  called  "ACTDUR"  calculates  the  hourly 
cost  at  optimal  machine  life  and  hourly  productivity  for  each  activity. 
These  values  are  used  to  calculate  activity  durations.  All  three 
values  are  outputted  with  a  description  of  each  activity  for  each 
set  of  standard  variables.  The  cycle  time  equations  contained  in 
subroutine  "PROD"  are  for  the  network  used  in  this  work. 

The  final  program  is  "GSPNA" .  The  programmed  versions  of 
Yen’s  and  Elmaghraby’s  algorithms  are  contained  in  this  routine. 
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C.2  Non-Standard  BASIC  Statements 

The  BASIC  programming  language  implemented  on  the  Hewlett- 
Packard  9845  desktop  computer  is  an  enhanced  version  of  standard  ANSII 
BASIC.  To  assist  readers  unfamiliar  with  H.-P.  BASIC,  a  list  of  non¬ 
standard  syntax  used  in  LOGNAP  and  an  explanation  of  its  purposes 
follows . 

Operators 


A 

Exponentiate 

& 

String  concatenation 

Statements 

COM 

-  Dimensions  and  allocates  memory  space  for  the 
specified  variables  in  a  ’'common”  memory  area. 

OPTION  BASE  1 

-  Specifies  that  the  lower  bound  of  arrays  begin 
with  one  rather  than  zero.  The  first  element 
in  array  A  is  A(l) . 

OVERLAP 

-  Sets  the  computer  to  an  overlapped  processing 
mode  to  increase  I/O  speed. 

REM 

-  Inserts  non-executable  remarks  (comments)  into 
a  program.  An  exclamation  mark  can  also  be  used 
as  a  comment  delimiter. 

WAIT  n 

-  Causes  program  execution  to  wait  the  specified 
number  of  milliseconds  before  it  continues. 

Functions 

INT (expression) 

-  Returns  the  Integer  portion  of  the.  evaluated 
expression. 

LEN (string) 

-  Returns  the  length  of  the  string. 

LIN (n) 

-  Outputs  a  carriage  return  and  the  specified 
number  of  linefeeds. 

PAGE 

-  Outputs  a  formfeed. 

RPT$ (string, n) 

-  Causes  the  string  to  be  repeated  the  specified 
number  of  times. 

SPA(n) 

-  Output  the  specified  number  of  blanc  spaces. 

TAB  (m) 

-  Output  the  next  item  starting  in  the  specified 
position. 

TRIM? (string) 

-  Delete  any  leading  or  trailing  blanks  from  the 
string . 

UPC$ (string) 

-  Returns  a  string  of  all  uppercase  characters. 
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Array  Operations 

MAT  matrix  =  matrix 

MAT  matrix  =  matrix 
(expression) 

MAT  matrix  =  ZER 

Mass  Storage  Operations 
ASSIGN  file  TO  //n 
ASSIGN  #n  TO  * 

CAT  TO  string  array 

CREATE  file,n,m 

LOAD  file 

READ  #1 , n; variables 


Copies  the  value  of  each  matrix  element 
into  a  second  matrix. 

Performs  an  arithmetic  operation  on  each 
element  of  a  matrix  with  the  result  becoming 
the  value  of  the  corresponding  element  in 
the  result  matrix. 

Sets  all  matrix  elements  to  zero. 


Open  a  file  for  accessing. 

Close  a  file. 

Copies  the  mass  storage  device  catalogue 
to  the  specific  string  array. 

Establish  a  data  file  of  n  records  of  m  bytes 
in  length. 

Place  a  program  stored  in  compiled  form  into 
memory. 

Positions  read  head  to  specified  record 
number  to  retrieve  variable  values. 


w., 

' 
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APPENDIX  D 

NETWORK  DEFINITION  AND  COSTING  INPUTS 


The  computer  output  contained  in  this  appendix  summarizes 
data  for  the  network  shown  in  Figure  7  that  was  analyzed  for  this 
research.  The  report  was  produced  by  the  INPUT  module  of  the  LOGNAP 
programs. 

The  first  section  of  the  output  summarizes  parameters 
which  in  effect  identify  the  network’s  size.  The  second  section 
outlines  the  linkages  within  the  network  in  terms  of  each  activity's 
initial  and  final  nodes  and  a  description.  The  final  section  contains 
the  data  specific  to  each  activity  which  was  used  in  the  calculation 
of  activity  durations. 
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APPENDIX  E 


COMPUTER  RUN  RESULTS 


Computer  output  summaries  of  runs  made  using  the  LOGNAP 
computer  programs  on  the  network  summarized  in  Appendix  D  follow. 
Individual  systems  of  equipment  were  ranked  in  ascending  order  of 
total  cost  per  cubic  meter  of  tree-length  wood  unloaded  at  the  mill. 
The  appendix  is  sectioned  as  follows: 

E.l.  Primary  Computer  Run 
E.2.  Sensitivity  Analysis  -  Fuel  Price 
E.3.  Sensitivity  Analysis  -  Average  Tree  Size 
and  Stand  Volume. 

Note  that  a  complete  computer  output  of  50  shortest  paths  is  included 
for  the  primary  computer  run.  In  the  sections  related  to  sensitivity 
analysis,  only  summaries  of  the  first  25  least-cost  paths  are  included. 

The  node  numbers  defining  the  shortest  paths  refer  to 
those  shown  in  Figure  7 . 

E.l.  Primary  Computer  Run 

Results  of  the  main  computer  run  are  contained  in  this 
section.  Note  that  these  results  were  also  used  as  the  standard 


conditions  in  all  sensitivity  analysis. 
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E.2.  Sensitivity  Analysis  -  Energy  Prices 

Results  of  the  computer  runs  made  to  determine  the  sensitivity 
of  forest  harvesting  systems  to  energy  price  variations  are  summarized 
in  this  section.  Table  E.l.  shows  the  fuel  price  assumed  for  each  run. 
Other  standard  variables  were  held  constant  at  values  used  in  the 
primary  run. 


Table  E.l.  Fuel  Prices  Assumed  in  the  Analysis  of  System  Selection 
Sensitivity  to  Energy  Prices. 


Fuel  Price 

($/L) 


%  of  Standard 
Price 


Computer  Run 
Number 


0.165 


75% 


RUN  1 


0.22 


0.775 


0.33 

0.385 

0.44 


100% 

PRIMARY  RUN 

125% 

RUN  2 

150% 

RUN  3 

175% 

RUN  4 

200% 

RUN  5 

1  ' 
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E.3.  Sensitivity  Analysis  -  Average  Tree  Size  and  Stand  Volume 

Results  of  the  computer  runs  made  to  determine  the  sensitivity 
of  forest  harvesting  systems  to  variations  in  average  tree  size  and 
stand  volume  are  summarized  in  this  section.  Table  E.2,  shows  the 
variable  values  assumed  for  each  run.  Other  standard  variables  were 
held  constant  at  values  used  in  the  primary  run. 

Table  E.2.  Average  Tree  Sizes  and  Stand  Volumes  Assumed  in  the 
Analysis  of  Optimal  Operating  Zones. 


TREE 

DBH 

(cm) 

STAND  VOLUME,  (m3/ha) 

60 

100 

140 

180 

12 

RUN  1 

RUN  2 

RUN  3 

RUN  4 

18 

RUN  5 

PRIMARY 

RUN  6 

RUN  7 

24 

RUN  8 

RUN  9 

RUN  10 

RUN  11 

30 

RUN  12 

RUN  13 

RUN  14 

RUN  15 
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